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Abo TRACT 


The  rediet  ion  from  sloe  antennas  on  a  cone  in  the  presence  of  an  inh  *»• 
gene ous  sheath  is  treated.  The  sheath  is  considered  cs  being  made  up  of  one 
or  two  conical  layers,  each  of  which  is  homogeneous.  The  boundary  conditions 
lead  to  a  system  of  integral  equations,  which  number  4M+4  for  a  sheath  composed 
of  M  (“  1  or  2)  conical  layers.  These  are  reduced  to  singular  integral 
equations  of  Cauchy  type,  which  are  solved  in  Iterative  fashion.  For  suffi- 
ciently  fine  stratification  of  the  sheath,  the  first  iteration  should  suffice. 

In  general,  fields  of  both  magnetic  and  electric  types  are  generated  in 
the  presence  of  a  sheath,  even  though  only  a  field  of  magnetic  type  may  be 
generated  in  free  space.  Foi  a  ring  slot,  however,  in  which  the  excitation  is 
axisaithally  symmetrical,  only  a  field  of  sugnetic  type  is  generated  even  in 
the  presence  of  a  sheath.  It  is  shown  that  the  solution  for  this  case  forms 
the  basis  of  the  solution  for  the  general  case. 

In  general,  evaluation  of  the  integrals  must  be  accomplished  by  contour 
integration,  which  leads  to  series  expansions  that  are  not  convenient  for 
numerical  evaluation.  For  thin  layers,  however,  Taylor's  series  expansions 
allow  all  but  one  of  the  coefficients  to  be  evaluated  in  closed  form. 

The  far  field  is  found  by  a  multi -dimensional  saddle  point  evaluation. 

This  is  illustrated  in  detail  for  the  free-space  case,  and  then  the  far  field 
patterns  in  the  presence  of  a  sheath  are  determined.  This  can  be  carried  out 
successfully  for  all  components,  and  to  arbitrary  orders  of  iteration. 

The  calculation  of  input  admittance  and  mutual  coupling  between  trans¬ 
mitting  and  receiving  slots  on  the  cone  is  formulated  and  methods  of 
calculation  are  discussed. 
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SECTION  1 


INTRODUCTION 

The  ionized  sheath  which  is  foraed  areund  a  high-speed  vehicle  or 
aissile  in  the  upper  ataosphere  has  a  profound  effect  on  the  transaission 
of  electrooagnetlc  waves  through  it.  This  is  aanifested,  for  excaple.  by 
the  blackout  of  transaissions  froa  a  space  vehicle  which  occurs  upon  re-entry 
into  the  ataosphere.  This  blackout  is  accompanied  by  a  large  change  in 
iapedance  of  the  on-board  transaitting  antenna  due  to  the  reaction  of  the 
aablent  aediua,  as  shown  by  in-flight  VSWR  aeasureaents .  The  antennas  on 
such  vehicles  are  typically  waveguide  arrays  of  resonant  longitudinal  slots 
The  vehicle  is  typically  conical  in  shape,  and  transaitting  and  receiving 
antennas  are  positioned  alternately  around  the  cone.  At  high  altitude  the 
sheath  is  conical  in  shape,  due  to  the  conical  torn  of  the  shock  wave. 

The  effects  of  a  plassa  sheath  on  antenna  radiation  and  propagation 
have  been  studied  intensively  for  slotted  cylinder  antenna*  [kef.  l]  and  to 
a  certain  extent  also  for  spheres  [Ref.  2}.  In  those  problems,  the  fonsula- 
tion  is  not  difficult  because  the  discontinuity  in  the  external  aediua 
occurs  in  the  radial  coordinate,  and  this  does  not  affect  the  separability 
of  the  wave  equation  or  orthogonality  of  the  functions  involved.  For  the 
conical  geooetry,  however,  orthogonality  does  not  exist,  so  that  the 
problea  be cooes  aore  complicated.  In  fact,  a  solution  of  the  electrooagnetlc 
boundary  value  problea  of  a  conducting  cone  covered  by  a  dielectric  or 
partially  conducting  conical  sheath  has  cot  appeared  heretofore.  This  has 
led  to  the  approxinate  representation  of  the  sheath  as  an  Infinitely  thin 
conical  impedance  discontinuity.  Pridaore-Brown  [Ref. 3)  treated  the  case 
of  a  magnetic  dipole  antenna  (loop),  end  Banos ,  et  al  [Ref.  4]  treated  the 
case  of  an  axial  electric  dipole  in  the  presence  of  such  a  sheath.  Their 
analyses,  however,  were  possible  only  because  of  the  infinitely  thin 
idealization  of  the  sheath,  and  for  Special  variations  of  sheath  impedance 
with  the  radial  coordinate  R. 

A  preoising  new  approach  to  the  solution  of  the  sheath  problea  was 
developed  in  a  previous  report  [Ref.  5].  A  net hod  of  using  the  Kontorovicn- 
Lebedev  (K-L)  transform  [Ref.  6]  was  proposed  to  solve  the  integral  equations 
which  represent  the  formulation  of  the  boundary  value  problea.  This  oechod 
is  developed  successfully  here,  leading  to  a  technique  for  solving  the  basic 
problea. 

In  Sec.  II,  the  basic  foraulation  of  the  problea  is  developed.  Although 
the  work  stateaent  of  the  present  contract  encoopasses  the  investigation  of 
single-  and  double -layered  sheaths  (two-  and  three -aediua  environments),  the 
foraulation  of  the  basic  equations  is  carried  ou t  for  a  sheath  c capos ed  of  K 
layers.  Consequently,  the  geooetry  considered  is  that  of  an  infinite  perfectly 
conducting  cone  provided  with  an  infiniteslsal  radial  slot,  the  cone  being 
overlaid  by  a  sheath  consisting  of  a  succession  of  conical  layers,  each  of 
which  has  an  arbitrary  complex  dielectric  constant  which  is  constant  throughout 
a  given  layer.  The  angular  thickness  of  a  given  layer  is  arbitrary,  as  is  the 
nuober  of  layers.  The  fields  in  each  layer  are  expressed  in  terms  of  sagnetlc 
and  electric  Hertz  vectors,  and  an  integral  representation  is  eaployed  for  the 
Hertz  vector  of  each  type.  Application  of  the  boundary  conditions  at  the  cone 
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and  at  each  layer  boundary,  together  with  the  requirement  that  the  field  in 
the  ambient  medium  be  finite  along  the  cone  axis,  then  results  in  a  system  of 
4(H+1)  coupled  Integral  equations,  where  M  is  the  number  of  layers  in  the 
sheath.  This  system  of  equations  constitutes  the  formulation  of  the  problem. 

In  Sec.  Ill,  the  simplest  problem  of  a  single -layered  sheath  is  considered. 

A  ring  source,  with  an  excitation  which  is  uniform  in  a2isuth,  is  considered 
first,  since  it  turns  out  that  the  solution  for  this  case  is  the  basis  of  the. 
general  solution.  this  type  of  source,  only  a  field  of  magnetic  type  is 

generated.  Each  integral  equation  is  converted,  by  means  of  the  K-L  transform, 
to  a  singular  integral  equation  of  Cauchy  type.  The  system  of  these  equations 
then  can  be  reduced  to  a  single  equation  of  Fredholm  type,  which  can  be  solved 
in  the  usual  manner  by  iteration.  It  is  shown  that  the  n£h  interacion  is  0(6n), 
where 

6  -  1-p3  ,  0  *  Ya  /Yi  • 

where  Yi  and  y2  are  the  complex  propagation  constants  in  the  sheath  and  in  the 
surrounding  medium,  respectively.  In  general,  the  evaluation  of  the  integrals 
in  the  solution  has  to  be  accomplished  by  contour  integration.  This  leads  to  a 
series  type  of  solution  in  terms  of  the  residues  at  singularities  of  the 
integrand.  These  include  the  zeros  of  the  Legendre  functions  as  a  function  of 
degree.  Accordingly,  e  coaputer  program  was  developed  for  Che  computation  of 
these  zeros.  The  analysis  simplifies  considerably  when  the  sheath  is  thin. 

For  this  case,  the  angular  functions  are  expanded  is  Taylor's  series  in  the 
angular  thickness,  of  the  sheath.  It  turns  out  that  the  first  approximation 
Is  G(^),  and  that  the  integral  equation  can  be  solved  In  closed  form  to  this 
order . 

The  case  of  a  sloe  source  is  taken  up  next.  This  involves  fields  of  both 
electric  and  magnetic  types.  The  system  of  Cauchy  type  integral  equations  is 
reduced  to  two  integral  equations,  which  may  be  called  the  excitation  and 
coupling  equations,  respectively.  The  coupling  equation  gives  the  coupling  of 
the  electric-type  field  to  that  of  magnetic  type  at  the  sheath  boundary,  while 
Che  other  expresses  the  exc 1 tat  Ion  of  the  magnetic-type  field  by  the  source, 
including  the  effect  thereon  of  the  electric-type  field.  I:i  the  case  of  a 
uniform  medium  (i.e.,  no  sheath)  the  electric-type  field  is  zero,  so  that  only 
s  magnetic -type  field  is  generated.  The  excitation  equation  again  may  be 
reduced  to  en  equation  of  Fredholm  type,  the  zero-order  term  being  Just  Che 
zero-order  term  for  the  ring  source.  From  this  zero  order  term,  the  electric- 
type  coefficient  can  be  determined  to  first  order,  which  is  0(6).  This  first- 
order  electric  coefficient  then  allows  the  reaction  on  the  excitation  of  the 
magnetic  coefficient  to  be  found,  which  is  0(6* ).  For  the  case  of  a  thin 
sheath,  the  magnet ic  coefficient  is  determinable  in  closed  form,  but  the  electric 
coefficient  Involves  an  integral  which  cannot  be  evaluated  in  closed  form, 
although  an  asymptotic  evaluation  is  obtained  later  in  Sec.  V,  in  connection 
with  the  determination  of  the  far  field. 

In  Sec.  J.V ,  the  extension  of  Che  analysis  to  the  case  of  a  two- layer 
sheath  is  worked  out.  As  in  Sec.  Ill,  this  is  carried  out  first  for  a  ring 
source.  Again,  an  exact  evaluation  of  the  resulting  integrals  requires  contour 
integration,  which  leads  to  a  multiple  series  type  of  solution.  For  (angularly) 
Chin  layers,  however,  the  integrals  again  can  be  evaluated  in  closed  form  for 
the  magnetic -type  coefficients,  but  not  for  the  electric -type  coefficients. 
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In  Sec.  V,  the  far  fields  and  patterns  are  calculated.  The  nethod 
eaployed  is  the  multi -disens ional  saddle  point  method  developed  by  Van  der  Pol 
and  Sremmer.  This  is  worked  cut  in  detail  for  Che  froe*space  case.  The  case 
of  a  single-layered  sheath  is  then  taken  up.  In  this  case  a  complication 
arises  because  the  integrand  possesses  poles,  a  situation  which  does  not  appear 
to  have  been  creaced  adequately  before.  A  method  for  dealing  with  this  situation 
is  developed,  and  applied  to  determine  the  fields  of  both  a  ring  source  and  a 
slot  source  in  the  presence  of  a  thin  sheath.  The  necessary  extension  for 
sources  of  finite  extent  is  also  worked  out.  Although  not  carried  out  in  detail, 
ic  is  pointed  out  that  the  technique  is  applicable  to  sheaths  of  arbitrary 
chickaess ,  as  well  as  to  the  case  of  a  mu It i -layered  sheath. 

In  Sec.  VI,  various  possible  extensions  of  the  net hod  of  analysis  developed 
In  this  report  are  pointed  cut.  The  calculation  of  input  admittance  and 
coupling  between  adjacent  slots,  both  of  which  were  formulated  in  a  previous 
report  [Ref.  5]  involves  the  near  fields,  so  that  asymptotic  expansions  cannot 
be  used.  The  evaluation  of  the  integrals  <sy  contour  integration  can  be  carried 
out  along  the  lines  eaployed  in  Sec.  III.  Alternatively,  numerical  evaluation 
of  the  integrals  would  be  required. 

See.  VII  summarizes  the  work  accomplished.  Various  mathematical  develop¬ 
ments  and  the  cocputer  program  for  the  Legendre  function  zeros  are  presented 
in  Appendices  A-D. 


3 


SECTION  II 


F0EKuLA7IC»  0?  IKTE(3LAI-  EQUATIONS 
2.1  GENERAL  CONSIDERATIONS 

For  a  conical  geaaetry,  a  spherical  coordinate  systen  (R  ,3 ,9)  is 
appropriate.  The  reduced  wave  equation  then  csn  be  separated  only  if 
the  electrical  properties  of  the  aedlus  do  not  depend  on  the  angular 
variables  6.9.  Hence,  for  an  implied  time  factor  eiu)t, 

cot;  depend  at  most  on  the  radial  variable  R.  An  infinite  cone  covered 
with  a  conical  sheath,  in  which  it  depends  at  s*ost  on  R,  is  thus  a  separable 
problem  in  spherical  coordinates.  Thus,  in  order  to  achieve  a  separable 
formulation  for  the  practical  case  where  the  sheath  properties  vary  wit  8, 
it  is  necessary  to  represent  the  sheath  as  a  succession  of  conical  layers, 
in  each  of  which  k  does  not  depend  on  8-  In  principle,  it  is  possible  to 
approach  a  continuous  variation  of  k  with  3  to  any  desired  degree  of 
appropriation  by  employing  a  sufficient  number  cf  layers. 

As  the  first  step  in  solving  the  sheath  problem,  the  case  of  a  sheat. 
whose  electrical  properties  are  invariant  in  the  radial  direction  will  be 
considered.  The  case  of  an  infinite  cone  will  be  considered  in  this  report. 

The  formulation  of  the  equations  for  a  sheath  consisting  of  M  uniform 
conical  layers,  as  in  Fig.  1,  will  be  carried  out  in  this  Section.  The 
solution  of  these  equations  will  be  carried  out  in  Section  III  for  the  case 
of  a  single  layer  (M*l) ,  and  in  Section  IV  for  a  double  sheath  (H“2) . 

A  general  type  of  field  nay  be  expressed  as  a  superposition  of  electric 
(TH)  and  magnetic  (TE)  modes,  derivable  froa  respective  Hertz  vectors  RH*.  Rli*. 
If  k  is  independent  of  R,  the  electric  and  magnetic  fields  are  given  by 


E"  *  K"1  iurl  curllKS1] 

E*  «  curl  iRV") 


i'  =  kturiqjX*) 

*  cart  curl  (RTT**) 


RUe,  RIF  each  satisfy  the  differential  equation 

(2.2) 

id*  d’t.lid  ill  v  '  - 

where  the  inhomogeneous  source  term  £  is  different.  In  general,  for  the 
electric  and  magnetic  nodes.  Solutions  of  (2.2)  then  v  a  be  expressed  in 
terms  of  solutions  of  the  corresponding  homogeneous  equation.  Separation 
of  variables  by  setting 

Rff  *  R  T{f)UCR)V(f) 

then  leads  to  the  differential  aquations 


jjHB  r>»*° 


(2.3a) 

(2.3b) 
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(2.3c) 


where  s3  is  the  separation  parameter  which  couples  (2.3a)  and  (2.3b),  and 
n3  is  the  separation  parameter  which  couples  (2.3b)  and  (2.3c).  For  a 
geometry  unbounded  in  the  qrcoordinate ,  as  in  the  cases  considered  here, 
n  is  an  integer. 

In  general,  a  solution  of  (2.2)  has  to  be  built  up  from  a  superposition 
of  solutions  of  (2.3a-c)  corresponding  to  a  continuous  range  of  the  separation 
parameter  s2  ,  leading  to  an  integral  representation.  A  method  for  evaluating 
the  coefficients  in  such  integral  representations  was  proposed  in  reference  5. 
The  basic  machinery  for  accomplishing  this  stems  from  the  K-L  transform.  This 
method  will  be  extended  in  this  report  to  attain  a  solution  of  the  boundary 
value  problem. 


2.2  FORMULATION 

The  antennas  of  interest  in  the  present  problem  typically  are  longitudinal 
slots  or  waveguide,  arrays  of  resonant  longitudinal  slots.  The  radiation 
problem  of  such  an  array  can  be  found  quite  simply  once  the  pattern  due 
to  a  6-function  slot  is  known.  If  the  field  of  the  latter  ar  a  typical  point 
is  denoted  by  E&,  then  the  field  due  to  the  array  is 

E  =  UiJ (2-4) 

Ti  a, 

where  f(RQ,%)  is  the  distribution  of  applied  field  over  the  array,  and  Ri  ,  It, 
and  cpL  ,%  are  the  bounding  coordinates  of  the  array.  Consequently,  the 
problem  is  essentially  solved  once  Ej  has  been  found,  since  only  the  integra¬ 
tion  of  known  functions  remains.  Therefore  the  analysis  to  follow  will  be 
concerned  with  finding  the  field  due  to  an  elementary  longitudinal  slot  at 
(Rq,^),  for  which  the  applied  field  is  given  by 

E,(8,)  =  E,8(R-V  *(*-*,)  (2.3) 

For  reasons  which  will  become  apparent  below,  we  now  choose 

»*  s(v-V*)(Vt’«sV*-)k  (2.6) 

For  convenience  in  the  applicatie-  of  the  K-L  transform  (see  Appendix  A), 
the  two  independent  solutions  of  th.  spherical  Bessel  equation  (2.3a)  that 
will  be  employed  here  then  are 

where 

*-ik  (2.7c) 


(2.7a) 

(2.7b) 


Iv  is  the  Bessel  function  of  the  first  kind  of  imaginary  argument 
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and  Ky  is  Che  Macdonald  function 


The  product 


„  1  ,(yR)-r,frR> 

M«)“  f - ^ - 


in  which 

R<  «  R,  R>  “  Rg  for  R  <  Rg 
R<  *-Rg  ,  R>  ■  R  for  R  >  K, 

then  is  the  appropriate  solution  of  (2.3a)  to  represent  a  6-funccion  source 
at  R  *  Rg . 

For  the  two  independent  solutions  of  the  associated  Legendre  equation 
(2.3b),  the  associated  Legendre  functions 

JD  *  Pj*_  (cos#) 

(2.9) 

%  -  Py?x(-ee*« 

will  be  employed  because  they  are  even  functions  of  v  and  the  first  is 
finite  at  0  *  0.  Again,  this  choice  is  made  for  ready  application  of  the 
K-L  transform.  These  functions  are  represented  here  by  f,  $  for 
compactness . 

The  trigonometric  functions  sin  no,  cos  ncp  will  be  employed  as  the 
independent  solutions  of  (2.3c). 

Consequently,  the  basic  form  chosen  for  the  Hertz  vectors  is 

Rf  *  i*  I  (2-10> 

4*  C 

where  <2  and  #  are  functions  of  the  integration  parameter  v,  and 

s'.  r.=  0 

«>0 

The  contour  C  and  the  choice  of  sin  rep  or  cos  rep  are  dictated  by  the  boundary 
conditions.  As  will  become  apparent  when  the  boundary  conditions  for  E^,  and 
Hq,  at  the  cone  and  sheath  boundaries  are  expressed,  cos  rep  will  be  required 
for  RIP1,  and  sin  rep  for  RIle.  Consequently  we  adopt  the  formulation 


RTI‘~  *  *v(*i*<)*  dv 

Rn^  =  £  $  (4;f  t  £;q),  ««  ,inn*vdv 


(2.11) 

(2.12) 


where  the  subscript  i  refers  to  the  icn  layer  of  a  multilayered  sheath,  and 
(  )y  denotes  that  all  quantities  within  the  parentheses  are  functions  of  v. 


2.3  BOUNDARY  CONDITIONS 

At  the  sheath  boundaries  0  *  0t  ,  the  tangential  components  of  the  fields 
must  be  continuous  At  the  cone,  the  applied  azimuthal  electric  field,  E^, 
is  given  by  (2.5),  while  the  radial  electric  field,  ER,  vanishes.  To 
complete  thl  notation,  to  the  angular  functions  -f,  \  we  affix  subscripts 
denoting  the  0-boundary.  Thus,  for  example,  at  0  *  0O  we  write,  in  accordance 
with  (2.9), 

■f.  =  fC/t(“»0.) 

i.  =  P„V'ce‘**> 

Corresponding  to  the  boundary  condition  (2.5)  at  the  cone,  0  *  0O ,  we 
obtain  from  (2.1)  and  (2.11) 

*  W  1  ft*  ®i  **  tn  co*  nf  v  dv  »  E  ,  Sftt-lO  (2.13) 

where  primes  denote  derivates  with  respect  to  0;  i.e.,  for  example. 

Corresponding  to  the  vanishing  of  E^  at  90 ,  we  obtain  from  (2.1)  and  (2.12) 

/(Ctf.  +  AV)*  =  0  (2.14) 

c 

(The  cn  is  superfluous  here,  since  the  n  *  0  term  vanishes;  it  was  used  in 
(2.10)  for  symmetry.)  Multiplying  (2.13)  by  cos  imp,  (2.14)  by  sin  n*p,  and 
integrating  over  <p  between  -tt  and  rr  leads  to 

*  tt*)"4  It  *(*->0  (2-15) 

iwtt.RO  ■*»(*,*>)  Wv  =0  (2'16) 

for  each  value  o£  g.  (2.16)  is  of  the  same  form  as  (2.15),  but  with  Ej,  *  0. 

From  (A3a)  of  Appendix  A  we  have 

SO?-fO*  -gr^r  T iv WR)A,WE,)vdi»*-jjTqgSs-6r04^*^)'>,<lv  (2.17) 

If  (2.17)  is  inserted  on  the  right  side  of  (2.15),  the  following  deductions 

may  be  made: 


(2.9a) 


(a)  The  contour  C  may  be  identified  with  the  imaginary  axis,  or  a 
contour  reconcilable  thereto  (i.e.,  without  crossing  singularities- of  the 
integrand) . 

(b)  The  integrand  of  the  left-hand  integral  behaves  properly  at 
infinity  of  the  Imaginary  axis  to  Insure  boundedness.  In  particular, 

(c)  y1  may  be  considered  as  real. 

(d)  From  the  K-L  transform  property  of  the  6-function,  the  integrand 
is  an  even  function  of  v;  this  means  that  ,  and  Bi  are  even  functions  of  v, 
since  jf  and  ^  were  chosen  to  have  that  property.  Thus,  it  is  not  necessary 
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to  distinguish  between  R<  end  R>  in  the  redial  functions;  either 


or 


ney  be  used  interchangeably  in  (2.11),  (2.12)  end  succeeding  integrals.  To 
show  this,  consider 

A  '  7t(v)i*0<)4t,(xjvdy 

where  E(v)  is  an  even  function  of  v.  By  replacing  K*(y  R)  in  (2.7b)  by  its 
definition  (2.8),  this  becoaes 

Changing  the  sign  of  v  in  the  first  tera,  we  obtain 

JL  =} 7 *  j$E6»)*,W*#C0sdy 

(e)  By  equating  integrands  on  both  sides  of  (2.15)  (or,  equiva¬ 
lently,  by  taking  the  K-L  transfora),  we  obtain 

(«.*'■ ♦*.*)  *£.  (2.18) 

where 

R.iTE,  (2.19) 

Similarly,  froa  (2.15), 

«,*♦*,*)  .O  (2.20) 

where  -i,  and  are  even  functions  of  v. 

At  the  sheath  boundaries ,  the  boundary  conditions  require  that  Eg,  %, 
E^),  be  continuous.  Writing 

%  *  ««.  / 

these  lead  to  the  respective  equations 

l  l*'-*  C  i<kf:  ♦  ii),  **(*>>  **<*.-)>*  *  0 

~  ^a;€ f  %  C\  ~  t/  }v4v  «  o 


(2.21) 

(2.22) 

(2.23) 

(2.24) 
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(2.25) 


♦*[&*♦*.  ***  =  0 


where 


-  Ai\/sin  8^  (2.26) 

<2-27> 

For  che  ambient  medium  in  an  M-layered  sheath,  which  includes  the  cone  axis 
S  »  0,  finiteness  of  the  field  requires  that  the  coefficients  of  ^  vanish. 
Hence 


<W.  =  o 


(2.28) 


(2.29) 


The  set  of  equations  (2.18),  (2.20),  (2.22)-(2.25) ,  (2.28),  (2.29) 
represent  the  formulation  of  the  boundary  value  problem.  The  method  of 
solving  che  integral  equations  to  determine  che  spectral  densities  <X  ,  8  , 

C  ,  <8  will  be  carried  out  in  Section  III  for  the  case  of  a  single  layer, 
and  in  Section  IV  for  a  two- layered  sheath. 
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SECTION  III 


SINGLE -LAYERED  SHEATH 


3.1  BOUNDARY  EQUATIONS 

For  a  single -layered  sheath,  the  boundary  equations  of  Section  II 
reduce  to  the  following  six  equations,  where  (2.28)  and  (2.29)  already 
have  been  applied: 


-i*  (3-1) 

«.f,  -  0  (3.2) 

i  [  j;  t  £,  v  H  =  O  0-3) 

)(.«-<,)[  (CUf.t  »,*),,  =  O  (3.4) 

«  *•  T»  13.51 


'  [(a,f^elq;),ivUOA<(y,3  -  (Qt A  fvrfv  =  O 


Since  the  right  side  »f  (3.1)  is  a  constant  given  by  (2.19),  it  is 
apparent  that  Ctv  ®ust  contain  a  factor  1/*  and  B,  a  factor  i/v* 
Siailarly,  fton  (3.2)  it  follows  that  t,  and  A,  ausc  contain  the  factors 
1/fw  .1/%..  respectively.  For  later  convenience,  therefore,  the  notation 
is  revised  lightly  at  this  point  by  defining 


A.  =  «LrO,i,ly.) 

A,  *  ICLtkiAld 

8,  *  (e,iji)w4j(y3 

C,  3  (t) * 

C,  »  ( ft),  Jhity,) 

D.  *  ^,V^*vly,) 

I.  *  £.J^(y.) 

In  addition,  we  introduce  the  paraaeter 

f  *  v»/V,  *  *t/4.  *yx/9, 


(3.7) 


(3.8) 


and  denote  Cj  by  x,  and  yx  by  y.  J_C_  will  be  assuaed  provls lonal Iv  in  the 
solution  of  the  integral  equations  that  p  is  real.  The  justification  of 
this  will  come  lit  the  eventual  Integration  over  R,  which  is  to  be  performed. 
The  boundary  equations  (3.1) -(3.6)  then  become 
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(A.i-B,)  =  £. 

(c,  t  D,)  =  o 

J  (»X-A1  ♦  6,^,/ti)  a»(0  -  (A1.£1),i.y(fx)3vdy  =  O 

J  {*.[•£  (C,r/f.+ 

*  C( A-fC^t  R.^).  i, U)  - CAx^a  }  * *»*  0 

J  ft  -  ti/CxfO,  j-*(px) 

t  N.CtA.f/f'  CAxf^CfO]}  -  0 

Proa  (3.1a)  and  (3.2a),  It  is  evident  that  the  source  generates  only  a 
ssgnetic-typc  field  at  the  slot.  Coupling  between  oagnctic-  and  electric- 
type  fields  occurs  only  at  the  sheath  boundary  through  the  E0~  and 
coeponents,  as  expressed  in  equations  (3.5a)  and  (3.6a),  respectively. 

A  special  case  of  interest  is  that  of  a  ring  source,  in  which  the  applied 
field  is  circularly  syaaetrical  around  the  cone.  This  special  case  will  be 
considered  first,  as  it  peraits  the  solution  of  the  integral  equations  to  be 
developed  in  its  sioplest  fora.  It  alr.o  foras  the  basis  for  the  general  case. 


(3.1a) 

(3.2a) 

(3.3a) 

(3-4a) 

(3.5a) 


3.2  RING  SOURCE 

For  a  circularly  syaaetrical  ring  swirce,  (2.5)  for  the  applied  field 
becomes 


Ef(6J  *  E.  KR-O  (3-9) 

The  field  is  then  independent  of  the  9-coordinate,  so  that  only  the  n  *  0  tera 
in  the  representation  (2.10)  is  required.  The  terns  containing  the  factor 
Nj  in  (3.5)  and  (3.6)  thus  drop  out,  so  that  no  coupling  between  electric- 
type  and  oagnetlc-type  fields  takes  places  at  the  sheath  boundary. 

Consequently  only  a  nagnetic -type  field  is  transaitted  through  the  sheath, 
as  in  the  free-space  case. 

Tims,  for  the  ring  source,  boundary  equation  (3.5a)  reduces  to 

J  lU.  ♦  R, */£), - (A»ff),  !,(**)] v  dv  *  0  ( 3  - 10) 

while  (3.2a),  (3.3a),  and  (3.6a)  do  not  apply. 

3.2‘.  1  Solution  of  the  Integral  Equations 

The  technique  for  solving  the  Integral  equations  (3.4a)  and  (3.19) 
is  based  on  the  X-L  transfora,  as  shown  in  detail  in  Appendix  3.  Taking  the 
K-L  transform  of  (3.10)  with  respect  to  x ,  we  obtain 
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0.11) 


(a,  f'M  *  =  "h  -f  dv 

where  u  Is  a  typical  point  on  the  iaaginary  axis,  the  contour  is  parallel 
to  and  to  the  right  of  the  iaaginary  axis,  as  shown  in  Fig.  B1  of  Appendix  B, 
and  M(n,v;p)  is  defined  by 

W  r.';C  <>•«> 

The  coefficients  ca(v,p)  are  given  by  (53)  of  Appendix  B.  h(u,v;p)  has 
singularities  at  v  »  -2a.  It  is  an  even  function  of  ii,  but  not  of  v. 

Siailarly,  the  transfora  of  (3.4a)  yields 

(3.13) 

{*jV*nd  ncfcr  ^e  eliainated  between  (3.1a),  (3.11)  and  (3.13).  The 

resulting  equation  aay  be  written  in  the  fora 


(O*  =  ^  K(^,v)  <piv)  dv 

(3.14) 

(3.15) 

mm  *  c --ft; 

(3.16) 

«r  *  w„/w, 

(3.17) 

*w-/w, 

(3.18) 

w«s  (f.v- 

(3.19) 

wi= 

(3.20) 

(3.21) 

t;=  (v*-w/(h‘-'A) 

(3.22) 

H(n,v),  which  is  an  even  function  of  both  u  and  v,  has  been  written  in  a 
fora  such  that 


Then,  by  defining 


*  I 


(3.16a) 


M,(p,v)  *  -  M(p,|»)  ~  M(n,»>)-|  (3.23) 

(3.14)  aay  be  rearranged  into  the  fora 

&K  t  jM,((i,v)qW)X(rl,w,p}dv  (3.24) 

*•  c, 

In  virtue  of  (3.16a)  and  (3.23),  the  integrand  of  the  second  integral  in 
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(3.24)  Is  not  singular  at  v  m  11.  Hence  v*  aay  shift  Cj  to  the  left  to  Cs  ,  the 
icaginary  axis,  in  the  second  integral. 

If  we  denote  the  second  integral  in  (3.24)  by  l(u). 

I<pl  "  "ijj  **  (3.25) 

It 

(3.24)  cay  be  written  as 

[SX~  !<*>  *  Tji  (3.26) 

This  equation  has  the  fora  of  (B6)  of  Appendix  B.  Hence  the  inversion  of 
(3.26)  is  gives  by  (BIO)  of  Appendix  5: 

*  4l  [l <U  -I(v)]Xt,..vtfMv  °‘27) 

The  first  tens  of  the  integral  in  (3.27)  aay  be  evaluated  as  follows: 

By  using  the  definition  of  I.  given  in  the  last  equation  of  (3.7)  and  the 
fore  of  HGi,v;l/s)  given  in  (B12),  we  have 

Interchanging  the  order  of  integrations ,  this  aay  be  written  as 

But  froa  (A3a)  of  Appendix  A,  the  v-integral  is  a  6-function.  Her«ce  we  obtain 

4UCp*>  *  (^.28) 


so  that  (3.27)  reduces  to 

9<r'  *  f.w  -  ±  h(v)  Xl^;**)  iv  (3.29) 

(3.29)  stay  be  solved  by  an  iteration  process  siailar  to  that  employed 
J'j  the  solution  of  integral  equations  of  Fredhoba  type.  For  I(v)  we  have 
froa  (3.25) 


I{»)  =  ^  (v,i,iq{A,)X(V<i^0d7. 

where  for  cOi  )  we  use  (3.29).  Then  we  obtain 

«**»•  *‘r»  -  MAW [q.00- 5; f IW Xu«v.$ 

*  *  ^1*^1  t 


(3.30) 


where 


<?.'H  1  ^7 ?  I X(e.v,y,/)dyjH)lv,»  «,C VX(A2;f)4* 
c.  ca 

?i'r'  * 

*  ^ 


(3.31a) 

(3.31b) 
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It  say  be  shewn  chat  CjjGi)  is  0(6n),  where 


*  *  \'f 


(3.32) 


This  follcvs  free  the  fact  that  I(^)  is  not  singular  at  v  *  n;  thus  Che 
a  *  0  tera  of  M(u,v;p)  in  (3.12)  gives  no  contribution  as  we  defora  the  contour 
to  cross  v  *  v*.  However,  froa  (,33b)  of  Appendix  B,  the  coefficients  cn(v,p) 
all  are  0(5)  for  a  >  0,  while  c,,(v,p)  •  c0(v,l/p)  *  1. 

3.2.2  First  Iteration 


For  the  first  iteration,  we  need  to  evaluate  (3.31a).  On 

using  (3.28)  for  ^(X),  this  is 

|*<r .* .»•*£*. (3.33) 

Frcn  (3.23)  and  (3.16) 

KUM*  C -  1  (3.33a) 

o/  and  orj,  froa  their  definitions  (3.17),  (3.18),  appear  to  have  poles  at 
tfj  =  0.  Htvever,  froa  (3.21)  Che  zeros  of  ^  occur  at  v-*i  «*  n.  3ut  then 
froa  (3. 19) 


K,  =  ?.(»*  tJ  Kl-cnt,)  -  &(.<•««,)  =  O 

since  ?-j( *cos*|)  *  (-)nPn(cosS) .  Sinilarly,  Kj  *  0  at  v-*f  *  n.  Hence  crj*  and 
ot%  have  no  singularities „  so  that  M.  (v ,A)  has  no  singularities  in  the  v-plane. 

The  inner  integral  in  (3.3)  will  be  considered  first.  Vith  U(v,X;p) 
in  the  fora  (3.12),  che  incegrand  converges  along  the  isaginary  axis.  In 
order  to  be  able  to  defora  che  contour  to  infinity  however,  it  is  necessary 
to  expand  A^(py)  in  accordance  with  (2.8),  and  then  defora  the  tera  with 
aj(py)  so  the  right,  the  tera  with  i.,(py)  to  the  left.  Alternatively  the 
sign  of  X  can  be  changed  in  the  second  tera,  whereupon  the  integral  beccoes 

The  contour  nay  now  be  closed  by  an  infinite  sea icircle  to  the  right.  Poles 
of  the  integrand  occur  at  the  zeros  of  (jtfj  in  K(v,X)  (see  (3.33a)),  at 
X  *  sv+2s  in  h(v,-X;p),  and  at  X  *  n.  Residues  at  X  *  n,  however,  can  be 
sh«*n  to  cancel  out  between  H(v,X;p)  and  H(v,-X;p).  Thus  we  obtain 


where 


d,  <») 


»l(V,f)  = 


v  -'Irj 


15 


(3.33)  then  becomes 


r  ^  J  V([».  dv 

A  contour  integration  in  the  v-plane  encounters  the  poles  of  K(p,v;l/o) 
of  The  result  is 

t  »»4  (  *•' r  *• 


f - *  . 

/  • 

•  \ 

*tv-# 

f  1 

>  \ 

V2 i 

l  « 

i  \ 

1  V2..n  h 
\[  1 

<  _\ 

(  »  * 

!  VI 

(3.34) 


where  >^(ar,3)  and  h(u.o)  have  been  def  ned  above. 


Since  Cj(p)  is  0(6* ).  an  iteration  stopping  at  c<  (p)  need  retain  only 
0(6)  terns  in  sp  (p)  •  Only  the  ^  »  0  tern  is  of  this  order.  Then,  caking 
use  of  (B3b)  and  the  expansion 


f  1  1  -  ♦  0(6*) 


we  obtain 


<f.(p5=  7  ^  rr~ 

%  r  5  ,  - ,  / _ i _ 1 _ .  1  !  \ 

L  S*^***'  V2*-*  A**2**?- 

♦  (j-$^  *  Vi^r)  *■  *  aS'} 

3.2.3  Second  Iteration 


(3-35) 


The  second  iteration,  c«(p),  is  given  by  (3.31b),  in  »fhich  c*  0-)  Is 
given  by  (3.34).  The  evaluation  of  Cg(p)  by  contour  integration  can  be  carried 
out  ip  a  manner  similar  to  that  used  in  the  evaluation  of  c^Ga)-  Oi  (p)  already 
contains  a  triple  swaaation  over  the  indices  k,a,p,  so  that  Ca(p)  involves 
three  additional  suxsacions,  caking  six  in  all.  It  is  evident  that  the 
successive  iterations  pose  a  severe  cooputational  problem  which  is  not 
attractive.  Thus  it  appears  that  it  nay  be  nore  economical  rroa  the  computa¬ 
tional  standpoint  to  divide  the  sheath  into  a  succession  of  layers  of  small  6 
to  Che  point  where  the  first  iteration  is  a  sufficient  degree  of  approximation. 

3.2.4  Thin  Sheath 

For  the  case  of  a  thin  sheath,  by  which  is  meant  one  for  which  the 
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quantity 


*.  *  A, 


Is  soall,  the  angular  functions  can  be  expanded  In  Taylor's  series  around  5  m 
It  turns  out  chat  the  first  approxioat ion  is  0(iJj  ) ,  and  that  the  integral 
equation  can  be  evaluated  in  closed  fora  to  this  order. 

Froa  the  Taylor  expansions  to  0 (d?  ) 

f,z  f.  '  >  f.'  »  r  >V  =  L 

f  -  {>  *  1 '  -  •>,  •»  y  d1  [coi’f,  ♦  c sc'd,  -  •/«,]  ■  f,' 


ve  obtain 


vhere 


Then  we  find 


*'» -  C  ‘  ~  2  *,>f*i*/*u*4 

Wt  /w,  -  i «  5,/va  5,  *  ■.  -  -3;  cct  f,  -  4  ’>* 


-  t  5:,V 

=  '£  $***-»»*)  ♦  O(tf) 


Then  (3.31a)  for  the  first  iteration  ci  (u)  becomes 


-  -i->1  )  d*  j  ^  ■*>  ifyi  ¥  ;*  <S 


(3-36) 


(3.37) 


Denoting  Che  inner  integral  by  J,  (v)  ,  as  before,  w*  write 

(>*-¥»)  - 

and  split  into  two  parts ,  to  be  denoted  by  a.  !>)  and  £-:»)  ,  respectively. 
Ford^alwe  obtain  directly  by  using  (BII)  of  Appendix  3  for  Ji(v,).;c)  and 
(A3e)  of  Appendix  A 


*  Nl-7*v  ~  J  =  (V1-  fJA/.y] 

For  d*U’)  ,  by  using  the  differential  equation  for  4j(py) 

*-  (*-*)*,(?*} 


(3.3S) 


(3.39) 


we  have,  siaiLarly, 


(3.40) 


* 


=  iKjp  -fjKty 
=  »*{§-!+  '-f*)^W 

on  using  the  differential  equation  for  fc*(y) ,  which  is  simply  (3.39)  with 
p  =  1.  Hence,  from  (3.38)  and  (3.40),  we  have 

i 

■i,N)  3  ^,ii^  -«£.!•; -  8yl.f<,ly)  (3.41) 

Inserting  this  for  the  inner  integral  in  (3.37),  and  using  (B12)  of  Appendix  B 
for  H(p,v;l/p),  we  obtain 

=  -■jAI&S«y2'Mpy)  (3.42) 

The  successive  iterations  now  could  be  written  dcwn  by  inspection.  But  the 
next  iteration,  <%(p),  in  virtue  of  (3.31b)  and  (3.36),  is  0(i3J ) ,  while,  from 
(3.36),  Oj  (p)  contains  an  additional  term  which  is  0(j5§),  and  thus  is  of 
higher  order  than  ^(p).  Thus,  only  the  first  iteration  can  be  retained  if 
Taylor  expansions  of  the  angular  functions  are  carried  only  as  far  as  O(^). 

3.2.5  Zeros  of  the  Legendre  Functions 

The  poles,  Xp,  which  arise  in  the  contour  integration  are  the  zeros 
of  the  Legendre  function 

(fo\  = 

More  generally,  in  the  case  of  an  unsymmetrical  excitation,  poles  occur  at 
zeros  of  the  associated  functions  Pvvl(c°£'90)  and  p£\£(cos90).  Accordingly, 
computer  programs  were  written  for  the  calculation  of  these  zeros.  These 
programs  are  described  in  detail  in  Appendix  D.  The  quantities  NU  and  MU  which 
appear  in  the  printout  are  defined  by 

*  0 

Consequently  the  poles  which  are  required  are  given,  respectively,  by 

XT  =■  NJ  -*  Vi 
X?  =  MU  +  ’/i 

3.3  SLOT  SOURCE 

i  We  now  return  to  consideration  of  an  elementary  longitudinal  slot  source, 

;  for  which  the  applied  field  is  given  by  (2.5).  The  formulation  of  the  integral 
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equations  has  been  given  in  Section  3.1,  where  it  was  pointed  out  that  coupling 
between  magnetic-  and  electric-type  fields  takes  place  at  the  sheath  boundary 
in  this  case. 


The  boundary  equations  now  are  (3. la)-(3.6a) .  Of  these,  only  (3.5a)  and 
(3.6a)  differ  from  the  ring  source  equations  considered  in  Section  3.2.  Since 
these  equations  contain  derivatives  of  the  radial  functions,  the  K-L  transform 
of  the  derivative  is  required  in  order  to  solve  these  equations  by  the  K-L 
transform  technique.  This  is  worked  out  in  Appendix  C.  The  applicability  of 
those  results  depends  on  the  condition  that  the  spectral  densities  have  no 
singularities  in  the  strip 

-(I4«l  <  Ra  v  <  |  +  « 

The  spectral  densities  turn  out  to  have  factors  1/.*,  and  1/flJ  ,  and  the 
derivative  occurs  only  for  angular  functions  of  order  n  >  0.  Since  the  only 
zero  of  these  functions  in  this  strip  occurs  for  ^  with  n  *  0,  the  condition 
is  met  in  the  present  problem. 

(3.1a)  and  (3.2a)  can  be  used  to  eliminate  B*  and  from  (3. 3a) -(3. 6a). 


The  latter  set  of  equations  then  becomes 

$  ■?*,(■*) (3.43a) 

c.  *» 

J  fJ(v).t1,i*)vdrf  -  i  dv  (3.43b) 

i  [7  -  fr(v)i>)]vdv) *  /tj (3.43c) 

C#  1  1 

i  *  i  E«itfw(v)i#(ps)  t  (3.43d) 

where 

W>  =  (d*-W  (C,  (3.44a) 

(■pC1)v  (3.44b) 

f,iv>  -  (V*-W  [(A,eS<)Ww/f;,;]/  (3 .44c) 

f4(v)=  (v*-'A)(f,AOv  (3.44d) 

ft(v)  -  N,  (C. (3.44e) 
f*(v)  =  M.  (f.C,),  (3.44f) 

f7<v)  *-[(*,  + (3.44g) 
f*M»  -(rfA,),  (3.44h) 

f,(v>  =  (3.441) 

=  (3.44J) 

f„(v)  =  'N,[(A,4S,)W^/ftV]v  (3.44k) 
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U.o  =  *.<.-**» 
u«'»  - 
Si  — 

5ts  M>.7w; 

Note  that 

=  inb>)/f4M  *  N/(*«-fe) 

The  transforms  of  (3.43a,b)  are 

c, 

respectively. 


(3.441) 

(3.44m) 

(3.44n) 

(3.44o) 

(3.44p) 


(3.45) 

(3.46a) 

(3.46b) 


From  (C6)  of  Appendix  C,  (3.43c)  and  (3.43d)  are  equivalent  to 

[Fs(v) *f7tv)]xJcx)  vdv  =  *'iy„(pAl *  ■f8Wlv(px)}vi’'  (3.47a) 

CiHf,(v)»F„(v)]Av(*)idy  =  ;  f (f^)iy(fx)  ♦  -ftIM C|*> + (3.47b) 
respectively,  where 


Fs,..W  *  *  (v'*>  (3.48a) 

**  =  (vI’A)/iy  (3.48b) 

The  transforms  of  (3.47a ,b)  are 

F,4t)  tf7(p)  =  fobM’**  *  ft(v)XCp,v;f)}dv  (3.49a) 

Ci** 

Vf>  t  Fi.(v>  *  rz  i  rfcMVfofcfl  +  X (3.49b) 

C*M 

where 

X*1*;  \  C  $}  X(p,v-l;f)  t  ^  XI ft#";')]  .  (3.50) 

In  virtue  of  (3.48a) ,  fru(pnj  are  required  in  (3.49a,b).  These  quantities 
can  be  obtained  In  the  following  way:  Multiply  (3.43a,b)  by  N,iK*,(x)/3? 
and  integrate  over  x  from  0  to  ®.  There  results.  In  virtue  of  (3.45), 

X(Ki|,^p)4v  (3.51a) 

Unt\)x-'k]  fI((KtO  =  f4WX(pitl,o;f)a>>  (3.51b) 


Since 
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we  obtain 


Hence 


where 


(«**-*' ^  *..*F*»>  =  it yt<nti,*;e)di 

r  <4*  T 

(K*-'/4)  F,.(^  -  g  I  Vv)  dv 

"A  e 


X**’  =  r't01'’A)X(p+>.‘>-,?)  +  (w*v*)W(fi-i,v;f)] 


(3.52a) 

(3.52b) 


(3.53) 


The  set  of  equations  (3.46a, b),  (3.49a,b),  and  (3.52a,b)  can  be  used  to 
eliminate  (Aj^ and  (C^  This  yields  the  integral  equations 

S  Xm-(h.v)  -  a  I  =<£),.  (3.54a) 


where 


C* 

(f.'AO, 
f'tvl  »  (fiO* 

«-w=£l£WG&$ 

"v<i  *  [’‘($U£),-  *(&(&]/(•-»  h  £l 

wrw» « n,(5^(^ 

*****  5  x^*  -  ^  *tM  =  I  »  ^rX(Kl,W;f) 


(3.54b) 


(3.55a) 

(3.55b) 

(3.55c) 

(3.55d) 

(3.55e) 

(3.55f) 


(3.55g) 


The  structure  of  (3.54a,b)  is  worthy  of  note.  (3.54b)  gives  c pe(v)  in  terms 
of  q/“(v) .  In  the  case  y^  “  ya  *  y  of  a  homogeneous  medium  (no  sheath) ,  p  *  1 
and  the  M6i±l,v±l; 1) ,  in  accordance  with  (B12)  of  Appendix  B  and  (A4a)  of 
Appendix  A,  become  6-functions: 

X{pU,v;l)  =  vji,  UjA^Wd*/**  =  Sfoif-v) 

Xtn.viljl)*  i(ft-v?l) 

so  chat 


Thus  for  Yi  “  Y2  •  ©e(n)  =  0-  For  Yi  FY2  »  Che  right-hand  side  of  (3.54b)  does  not 
vanish,  so  that  then  oe(u)  does  not  vanish.  Thus  (3.54b)  expresses  the  excitation 
of  the  electric-type  field  by  the  magnetic -type  field  at  the  sheath  boundary. 

In  (3.54a),  the  second  Integral  similarly  vanishes  for  Yi  *  Ya  »  30  t*iat 
the  equation  then  expresses  the  excitation  of  the  magnetic-type  field.  For 
Yii^Ya*  the  second  integral  represer-ts  the  alteration  in  the  excitation  of  the 
magnetic-type  field  due  to  the  creation  of  the  electric-type  field  at  the  sheacn 
boundary. 

3.3.1  Solution  of  the  Integral  Equations 

The  solution  of  the  integral  equations  (3.54a,b)  can  be  effected 
in  a  manner  similar  to  that  employed  in  Sec.  3.2.1  for  the  case  of  a  ring  sou rce. 

As  can  be  seen  from  (3.55c),  ,v)  in  (3.54a)  is  the  same  as  the  function 

K(p,v)  which  was  obtained  in  the  ring  source  case,  and  given  by  (3.16),  so  that 
lf*(n,li)  *  1.  Since  <p*(}i)  is  zero  for  6  *  0  (p  »  l),  it  folltxre  that  o®(v)  is 
at  most  0(6);  i.e.,  must  lead  to  terms  that  are  0(6)  at  most.  In  fact, 

it  is  shwn  in  Appendix  B  that  Jf5*1 leads  to  terms  which  are  0(6).  Then  it 
follows  that  the  second  integral  in  (3.54a)  must  be  0(63).  Thus,  to  0(6),  the 
generation  of  the  electric-type  field  at  the  sheath  boundary  does  not  affect 
the  excitation  of  the  magnetic-type  field. 

In  view  of  this  property,  we  can  define 

=  M,(pl(i)  =  0  (3.56) 

and  write  (3.54a)  in  the  form 

J-  Sq-{v)35(Hlv;pldv  =  (g^- l"{K>  t  J»(H)  (3.57) 

M 

where 
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(3.58a) 

(3.58b) 


rip)  =  ^  5  <TM  XCy.Vtf)  dv 

c» 

e,«t 

Ia(u) ,  as  in  Che  ring  source  case,  is  0(6).  As  pointed  out  above,  J*0O 
is  0(6*). 

As  in  Sec.  3.2.1,  Che  inversion  of  (3.57)  is 

f'lp)  *  J (3*59) 

1  c, 

(3.59)  can  be  evaluated  iteratively,  jusc  as  in  Sec.  3.2.1.  In  this  procedure, 
the  results  differ  froa  those  obtained  there  in  only  one  essential  respect, 
namely  the  contribution  to  cP(u)  froa  Jc(v) . 

In  (3.54b),  M*0i,v)  has  been  defined  in  such  a  way  that  ^(ji.p)  *  1. 

Hence  we  can  define 


(3.60) 

and  write  (3.54b) 

in  the  fora 

^  J  f'M  =  -  I *(y)  *  J» 

(3.61) 

where 

I V  =  i 

(3.62a) 

rQ»)  *  5  <(v) 

M* 

(3.62b) 

(3.61)  differs  in 
right-hand  side. 

fora  froa  (3.57)  only  in  the  absence  of  a 
The  inversion  of  (3.61)  thus  is 

source  tera  on  the 

4 

(3.63) 

The  "source"  of  oeOO  is  the  second  tera,  JB(v) ,  on  the  right,  which  steas 
from  cji*(v)  in  virtue  of  (3.62b).  Due  to  the  property  of  this  integral 

is  0(6).  Then  the  integral  of  the  first  tera,  Ie(v) ,  is  0(6*)  in  virtue  of 
(3.62a). 

If  we  express  cPOO  and  9*00  as  a  sua  of  iterations  of  successively  higher 
orders  of  6,  i.e.. 


Wit]  (3.64a) 

r  i*a  ‘ 

(3.64b) 

where 

•?;%>  =  olz>) 

then  (3.58a,b)  and  (3.62a,fc)  beccoe 
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i*V)  =  £  iJV -  i  a  lnTwif!,(piXiK*ifld*  (3.65*) 

J*<p)  =  t  3-V  *  £  -M  (3.65b) 

j*1  j*i  t,« 

I*(}d  *  2  Xfp,g;p)diJ  (3.65c) 

J*40  I  J;V  =  2.  ^  ImJ  W  X1*^  (3.65d) 

respectively.  As  pointed  out  earlier,  the  IBGi)  and  I*(u)  integrals  are  0(6) 
relative  to  the  ©-function  in  the  integrand  because  ^“(p.p)  m  0  m  Hf(p  ,p) . 

Thus  IB(p)  and  I*(p)  are  0(6^).  Sinilarly,  because  of  the  property  of  )ru,lv, 
J*(p)  and  JB(p)  are  0(6*).  Thus  it  is  evident  fro a  (3.59)  that  p®(p),  JX), 
arises  from  1®,  i.e.,  ©JLj  ,  and  J*.  But  froa  (3.65b)  J*  steas  froa  : 

Hence  cj’GO  cooes  froa  ccjL*  and  cp_i  . 

Siailarly,  froa  (3.63)  it  follows  that  ©*(p)  arises  froa  I*  and  J®,  i.e., 
frcn  of.\  and  p?_i. 

Since  p*(p)  has  no  source  tens, 

=  O 

Then  it  follows  that  q*f(;i)  and  cf(p)  are  given  by  exactly  the  saae  expressions 
as  for  the  ring  source  case  considered  in  Sec.  3.2.  The  only  difference  in 
these  terns  between  the  ring  source  and  the  slot  source  cases  is  that  for  the 
ring  source  the  azisuthal  parameter  n  of  (2.10)  is  0,  while  for  the  slot 
source  all  values  of  n  are  involved.  Thus  froa  (3.28)  we  obtain 

Wpl  =  (366> 

while  c£*(p)  is  given  by  (3.34),  (3.35),  or  (3.42),  or  in  integral  fora,  by 
(3.32). 

Correspondingly,  ©f(p)  is  given  in  integral  fora  by 

*  (3.67) 

The  evaluation  of  (3.67)  by  contour  integration  can  be  carried  out  as  in  the 
evaluation  of  (3.32).  In  (3.67),  takes  the  place  of  X(v,X;p)  in  (3.32), 

and  Mf(v,X)  appears  instead  of  ^(v.X).  The  poles  of  X(v,X;p)  which  gave  rise 
to  residues  at  X  *  isrr2a  now  are  supplanted  by  the  poles  of  at 

X  *  ±v+2a+l.  While  there  was  no  pole  in  (3.32)  at  a  *  0  because  K  (v,v)  ■  0, 
there  is  such  a  pole  in  (3.67)  because  Ha(v,v)  is  not  zero  for  p  f  1.  The 
poles,  Xp,  of  K  (v,X)  occurred  at  the  zeros  of  ,  while  the  poles  of  fcf(v,X) 
occur  at  the  zeros  of  as  can  be  seen  froa  (3.55d).  Except  for  these 
differences,  the  evaluation  of  ©f (p) ,  and  higher  iterations  if  desired,  can  be 
carried  as  in  Sec.  3.2. 

3.3.2  Thin  Sheath 

In  Sec.  3.2.4,  the  thin  sheath  approx inat ion  was  introduced  by 


defining 


end  developing  the  Legendre  functions  in  Taylor’s  series  around  @0 .  The  first 
iteration,  which  is  0(6),  then  turned  out  to  be  O(d^)  and  could  be  evaluated 
in  closed  fora.  The  second  iteration  then  would  be  0(t)f),  so  that,  if  the 
Taylor's  series  developaent  is  carried  nc  higher  than  iJj  ,  one  is  necessarily 
restricted  to  the  first  iteration.  A  siailar  procedure  will  be  followed  here. 
In  particular,  teras  will  be  United  to  the  lowest  order  in  . 


Proa  (3.36) 


Sinilarly,  we  find 


m;w)  fO(**) 


It  has  already  been  oointed  out  chat  ©e(\>)  is  0(6).  Consequently,  if  a  solution 
is  limited  to  0(6)  tarns,  it  is  evident  froa  (3.54b)  that  H*(|i,v)  must  be 
United  to  0(1)  ternr .  Hence  we  take 


so  that,  in  accordance  with  (3.60) 


MfCp.vl  «  fcSfU-K*)  * 


(3.68a) 


«rw>  =  Hi  ♦ 

M{(M  (fit 

Of  the  iterative  teras  in  (3.59)  for  q?*(n) ,  we  then  find 


IS(V)  -  0(3  6) 
Je(v)  -  0(6*) 

while  in  (3.63)  for  ^(ji) , 

I«(v)  -  0 ($6*) 
JB(v)  -  0(^  6) 


froet  (3.58a)  and  (3.68a) 
from  (3.58b)  and  (3.68b) 


froa  (3.62a)  and  (3.68a) 
froa  (3.62b)  and  (3.68c) 


As  in  Section  3.2.4,  we  then  find 


(3.68b) 

(3.68c) 


(3.69) 


qfGO  steas  froa  Ja(v)  in  (3.62b).  This  integral,  however,  cannot  be  evaluated 
in  closed  fora,  so  chat  we  have 
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(3.70) 

(3.71) 


-  ~~rT  S/.(r) 

C,»l 

<W *  S3- UV#M^ 

Although  (3.71)  cannot  b«  evaluated  in  closed  fora,  nevertheless  it  is 
possible  to  obtain  saddle  point  developnents  of  the  far  field  of  electric  type 
(i.e.,  stealing  froa  cpf).  This  will  be  carried  out  in  Sec.  5.4. 


26 


SECTION  IV 


DOUBLE-LAYERED  SHEATH 


4.1  INTRODUCTION 

The  formulation  of  Che  boundary  equations  for  a  sheath  composed  of  H 
uniform  conical  layers  was  given  in  Section  II,  and  in  Section  III  the  case 
of  a  single-layered  sheath  was  treated.  In  this  Section,  a  double -layered 
sheath  will  be  treated. 

For  an  M-layered  sheath,  the  two  source  equations  (2.18)  and  (2.20)  at 
the  cone  surface,  the  four  equations  (2.22)-(2.25)  expressing  the  continuity 
of  the  tangential  components  of  field  at  each  sheath  boundary,  plus  the  two 
equations  (2.28),  (2.29  expressing  the  finiceness  of  the  field  along  the 
cone  axis  in  the  aabien;  aediua,  sake  up  a  totality  of  4(H+1)  boundary 
equations.  After  introducing  the  notation 


\  m  °-i  JMy;) 

B.  JUtf 

c*  * 

(4.1) 

these  equations  become 

(A,^e  e,tA  5  c. 

(4.2a) 

*  O 

(4.2b) 

/  lv*-W[7(C4^Rsv)^  '  7-CCu.f.r  =• 

C  '*  •*! 

0  (4.2c) 

Ui*-/4  «•  4j,*0OJvd»-  =  0 

(4. 2d) 

(4-2e) 

5  IT; (C;*'* $.<<). r *<), V*.-..) 

<4.2f) 

The  solution  of  the  problea  of  a  single -layered  sheath  was  developed  in 
Section  III,  first  for  a  circularly  sysraecrical  ring  source,  and  then  for  an 
elementary  slot  source.  It  turned  out  chat  Che  solution  for  the  ring  source 
fores  the  basis  for  the  sore  general  case;  that  is,  the  ring  source  solution 

is  the  "zero  order"  solution  for  the  general  case.  For  the  culciple-iayered 
sheath,  Che  case  of  a  general  source  distribution  also  is  an  extension  of  the 
ring  source  solution.  Consequently  the  case  of  a  ring  source  will  be  considered 
first,  in  particular  for  a  double  sheath  (M*2). 
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4.2  RING  SOURCE 


For  a  circularly  syrsaetrical  source,  there  is  no  coupling  at  the  sheath 
boundaries  between  nagnetic-type  and  electric-type  fields.  Hence  for  a  ring 
slot  excited  by  an  azinuthally-directed  electric  field,  only  a  nagnetic-tj pe 
field  is  created.  Thus,  the  terns  in  (4.2e,f)  which  contain  derivatives  of 
the  radial  function  drop  out,  and  the  boundary  equations  (4.2b,c,f)  becoae 
superfluous.  (4.2e)  then  reduces  to 

=  I  (4.3) 

Follcwing  the  procedure  used  in  Sec.  3.2.1,  on  the  right  side  of 

(4. 2d)  and  (4.3)  is  expanded  into  a  series  in  iv.2*(xi)  an<*  c*le  K"L  transform 
of  the  equation  taken  with  respect  to  xt .  Froa  the  transforn  property  discussed 
in  Appendix  A,  each  x{  is  considered  as  real.  This  gives 


(m*-W  Wiffe  ♦  5  * 

(4.4a) 

c  51  $  *'*  K*’  ** 

c% 

(4.4b) 

where 

*r  * 

(4.5) 

Pi  =  I*. /rL 

(4.6) 

Again,  each  pt  is  considered  as  real,  since  ultinateiy  an  integration  over  the 
source  coordinate  (i.e.,  the  yt )  can  be  taken  along  the  contour  for  which  ot 
is  real. 

By  means  of  (4.2a),  B1  nay  be  eliminated  from  the  i  *  1  equation  of 
(4.4a, b),  and  then  Aj  may  be  eliminated  from  the  resulting  two  equations.  This 
leads  to 

(4.7) 

where  is  defined  by  (3.7),  and 

«.w 

(4.8a) 

K-V**  r  (Si  (?fi  (-£),  *  I 

(4.8b) 

The  remaining  equations  of  (4.4a,b)  (that  is,  for  i  *  2)  are  solved 
simultaneously  to  eliminate  (3^,..  The  result  can  be  expressed  in  the  form 

*xV  * * 5T 51 * «?(*/) <Jw J xr'dv 

(4-9) 

where 

’fpsaia-fta® 

(4.10a) 

(4.10b) 
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(4.20) 


*;(",V)  =  X  ( n  ;  Vp-) 

Since  lf,b(p)  are  0(6^,  if  we  write 


<?;v « 

JiO  »i 

(4.21a) 

*‘V  -  £ 1?*H» 

(4.21b) 

where 

<•114=01^  i;-‘(w=0(Si) 

(4.21c) 

then  = 

0.  so  that  from  (4.19) 

W*1'  =  dv 

(4.22a) 

In  particular, 

in  virtue  of  (4.13b), 

<?»,„(**>  r-L  J  ^(v)  =  ^U<(y)X(H.v;'/f,)dP«  Mh(M> 

(4.22b) 

Th'.s  result  may  be  inserted  in  (4.22a)  with  i  =  2  to  yield 


=  MmU 'tl) 

Hence 

<»..  Cf^  =  ;*'.2  (4.23) 

For  i  =  2,  it  follows  from  (2.28)  that  B3(v)  =  0,  so  that  c^(\>)  =  0. 

Hence  from  (4.16)  with  i  =  2  and  (4.13a) 

-V!  lp>®  &} s  -  <??«•►*!  (4.24a) 

from  (4.22b).  Obviously,  then,  q^,0(u)  =  0.  Then  from  (4.16)  with  i  =  1  we 
obtain 

—  f o lv* d>> -f.l0lj«)=  5 =  $£f<fd  (4.24b) 

so  that  cf0(u)  =  0.  Hence 

<?*„{«)  r  O  1,2,3  (4.25a) 

Thus  gj’fa)  is  0(6j).  (4.22a)  then  becomes 

'  44  1=1.2  (4.25b) 

(4.25a,b)  are  a  result  of  the  fact  that  lf,b(p)  =  0(6S)  relative  to  csf+,b(v) . 

The  generation  of  cpb(u)  takes  place  via  (4.13b),  which  can  be  written  as 

H,<jO  =  r  £(<??,;  (pH <?,*(*») ]  *  f.iuly) 

j  *0  * 

But  from  (4.24b),  £,-!,,(>’)  is  just  the  value  of  qfi0(p).  Hence  it  follows  with 
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(4.25a)  that 


so  that 


Vi(u)  =  ° 


I' '.2.  • 


Thus,  from  (4.19)  for  i  =  1 


1  57  /  C$;,(v)  -  I,1  (v)  -  X‘ (vi]  B^dv 


(4.26a) 


(4.26b) 


(4.27a) 


where  if’^(v)  are  given  by  (4.17).  Subtracting  out  the  j  *  0  terms  from  this 
equation  by  means  of  (4.25a,b'>  and  using  (4.17),  we  obtain 

£  fij  <t*l  =  -  £  J(  1,*M  r  ifMjX.^dv 

j.t  a 

*  -  7%  IS,'"'**  I fx*U)  *  (4 . 27b) 

«*  <=•  „  v 

<lyO 

The  j  =  1  term  is  obtained  by  using  0(X)  in  the  integral.  But,  since 
«2(X)>  by  (4.25a),  is  0(63),  the  I^(>.) ’integral  is  0(6|),  so  it  does  not 


contribute  to  the  j  =  1  term.  Then  (4.27b)  for  j  =  1  becomes 


S j W» (4.27c) 

where  is  the  portion  of  which  is  0(6S);  similarly  for  5^’^.  Froa 

(B7)  and  (B3a,b)  of  Appendix  B,  for  a  general  value  of  i, 


(4.28a) 

(4.28b) 


in  virtue  of  the  relations 


h*=  l 

1  -  Yf.'  =  -  i:  v  OU/1 

Next,  from  (4.19)  for  i  =  2,  we  get  in  a  similar  way 
"  itT  dv  = 

since  c£(p)  =  0  from  the  boundary  condition  (2.28). 


(4.29) 


Thus  all  the  Oi.,i(p)  have  been  determined.  However,  except  for  i  =  1 
and  i  =  3,  the  components  ©f 1 j  and  ©ftl  are  yet  to  be  found;  i.e.,  we  need  to 
find  For  these,  it  is  necessary  to  solve  (4.9)  and  (4.11)  separately. 

This  can  be  carried  out  in  exactly  the  same  way  that  was  used  in  the  solution 
of  (4.14).  Thus  we  define 


-  Mf(p,V)-l,  M;*(c,m)  =  0 


(4.30a) 


K ‘  (H,v>  =  =  Kj’W- 1,  *  0 


(4.30b) 


Note  that  we  already  have  from  (4.10b)  and  (4.12a) 


=  0 

For  notational  convenience,  we  can  write  subsequently 


Ma  (jj.v)  =  M.  (.»,>>) 


so  that 


-  Ki  (w.vi 


Then  (4.9)  and  (4.11)  may  be  written  as 


lr  i>  =  -  Jf(fl)  -  Jjfy 

i*Lo'*r'<»  *  fV  -  jfw  -g.*4o 


respectively,  where 


3?*W  =  ^  i  <?.?(»'  iO'dv 

gv* ?! !  *;%.•>> 


(4.31a) 

(4.31b) 


(4.32) 


(4.33a) 

(4.33b) 


(4.34a) 

(4.34b) 


In  virtue  of  (4.32),  the  J{  and  frj  integrals  are  non-singular,  so  that  they 
are  0(6t).  Hence  (4.33a,b)  may  be  inverted  to  yield 


(4.35a) 

(4.35b) 


-  ^lv>  -iktAtll"**  (4  •  35b) 

^  a  ,b 

The  zero-order  approximations,  j  =  0,  of  c.,j  have  already  been  found  in 
(4.25a,b).  Hence  the  j  =  1  approximations  can  be  evaluated  directly  by 
subtracting  out  the  j  =  0  terms  of  (4.35a,b)  to  yield 


lU,  ^  --  71  [  tCM  '  J;>>  *  *AM3Sfd* 

O'*  •  i  S  "9 1..W35;"'* 


(4.36a) 


(4.36b) 


a.b  a  ,b  ^  a,b 

where  Jt  x  and  j  are  given  by  (4.34a,b)  with  In  particular, 

for  i  =  2 


^  =  71  ! C  <v»  -  Jv  Ml 

<?>>  =  7T  S[?l«  -  &{v)  -  gi/vJjH-'dv 


(4.37a) 

(4.37b) 


Since  =  0  from  the  boundary  condition,  $j(v)  =  0  from  (4.34b),  so  from 
(4.37b)  and  (4.34b)  it  follows  that 


=  9v  V)  8  -57 


(4-3S) 
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Hence 


(4.39) 

But  already  has  been  evaluated  in  (4.29a)  for  i+1  =  9  a 

is  determined  by  (4.39).  2‘  Hence  ^.l^) 

Cl!e  "B*  °f  PEocedufe  can  be  followed  to  evaluate  the  j  =  2 
terms.  In  this  I?,2(n)  and  „m  first  appear,  since  h“  ™ J  2 

‘“r1''''™* hliu) 

~st  °f  the  is  - * 

4.2.1  Thin  Sheath 

i-  X-ylor'.  series  ^  «*  ^eI^d 


^=5.-6; 


(4.40) 


Furthermore,  this  approximation  is  limited  to  the  first  iteration  or  O ft  1 

(i%  ,^l"2  0<5>>.  ifto).  J-%)  •  !,b(u)  text*  Jeflbed^ 

„ce;''t“  «■*“■»,  respectively,  ,«  0(6?),  so  that  these  terns  can  be 

«-  “g-  *  <«•»>  -  O-UO.  ve 


o.*u, 

TO.'P.v)  =  (-^)r(^,  -r)  -I  (1.41b) 

where  Wt0  ,  W,0  ,  and  Ul  are  obvious  generalizations  of  (3  19W3  211 
analogous  to  (3.39),  we  find  '  *'  '3*Zl'*  Then> 

(4.42) 

(4.42)  does  not  involve  the  angular  functions,  so  that  all  the  integrals  can 
(4.25b)Ua^d  f rocjC (4 1 34a) ° for ^xamp le?*  ^  SeCti°n  3'2'*’  Fr°°  (4'1?)  3nd 

=  izll  (4.43) 

which  is  the  counterpart  of  (3.45). 

(4.23Ke,Slj4!mf"t1herar«°8^b'“l,“'1“'  Fr“ 

<m)  =  — j  SKt  ly.W  ■*»(?;?> (4.44) 
where  Kffc.v)  is  defined  by  (4.12a).  Due  to  the  form  of  Kffo.v),  this  integral 
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cannot  be  evaluated  in  closed  form.  However,  it  is  possible  to  obtain  a  saddle 
point  development  of  the  far  field  involving  this  integral,  in  the  manner 
developed  in  Sec.  5. 


SECTION  V 


FAR -FIELD  PATTERN 

5.1  INTRODUCTION 

In  this  Section,  Che  far-field  pattern  of  a  sheath-covered  slot  antenna 
will  be  determined.  The  method  eaplcyed  is  the  saddle  point  aethod,  due  to 
Debye,  which  is  frequently  used  for  the  asynptotic  evaluation  of  integrals. 

Van  der  Pol  and  Brenner  [Ref.  8]  extended  this  aethod  to  the  evaluation  of 
a  mu It i -dimensional  (i.e.,  multiple)  integral  encaintered  in  the  diffraction 
of  radio  waves  around  a  sphere.  The  technique  requires  that  all  functions 
which  occur  in  the  field  representation  be  expressed  as  exponential  integrals. 

A  Taylor  series  expansion  of  the  exponent  to  second  order  is  made  about  its 
stationary  point ,  or  saddle  point,*  whereupon  the  integral  becomes  a  multiple 
Fresnel  integral.  In  the  present  problem,  a  complication  is  encountered  because 
the  integrand  possesses  poles,  the  poles  being  those  which  occur  in  the  function 
K(u.v;p),  defined  in  (57)  of  Appendix  B.  This  situation  does  not  appear  to 
have  been  treated  in  the  literature  heretofore.  A  method  for  dealing  with  this 
situation  is  developed  here. 

As  an  illustration  of  the  aethod,  the  derivation  of  the  far  field  of  a  ring 
source  on  a  cone  in  free  space  by  the  multi-dimensional  saddle  point  aethod, 
which  does  not  seen  to  have  been  determined  by  this  method  previously,  will  be 
developed  in  detail  first.  This  method  will  then  be  extended  to  the  case  of 
a  sheath-covered  slot.  This  will  be  carried  out  first  for  a  ring  source,  in 
which  only  a  field  of  magnetic  type  is  set  up.  Following  this,  the  case  of  a 
slot  source,  where  fields  of  both  magnetic  and  electric  type  are  generated, 
will  be  worked  out.  The  far-field  patterns,  for  an  infinitesimal  slot,  as  well 
as  a  half-wave  slot,  will  then  be  determined.  The  extension  to  an  array  of 
slots  will  also  be  given. 


5.2  RING  SOURCE  IN  FREE  SPACE 

For  a  ring  slot  source  in  free  space,  excited  by  a  circularly  syaMetr'  al 
azimuthal  electric  field  Eg,  the  only  electric  field  component  is  E  .  Th-  _ 
from  (2.11),  (2.21),  and  (3.1)  with  n  -  0, 

Et  =  i  “  "S  |  (aA‘,W  v  dv 

im  . 

~  ~  iJAkW  v  rfv  (5.1) 

Introducing  the  defining  relations  (2.7a,b)  for  ji,,  JL,,  this  may  be  written  as 

-RE*  =  {RR,)a  £,*  (5.2a) 

where 

4  =  S  I, W  K,(y)»dv  (5.2b) 

*  The  terms  "saddle  point"  and  "stationary  point"  will  be  used  interchangably 
hereafter. 
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For  the  application  of  the  saddle  point  nethod,  integral  representations 
are  required  for  the  functions  in  the  integrand.  For  the  Iv  and  K*  we  have 
the  Sooner fe Id  integral  representations 

M*>  *  (5.3a) 

—  T  . (5.3b) 

For  the  angular  functions,  it  is  desirable  at  this  point  to  replace  the 
Legendre  function  c m  Pv-^(cos8)  by  the  traveling  wave  angular  functions  used 
by  Felsen  [Ref.  9  3  ,  which  are  defined  by 

where  Pv-^,  Qy-%  ar<  the  usual  Legendre  functions  of  zero  order.  Hence 


f  * 


have  the  Laplace  integral  representations 

yn,{U  _  -  1  y  [C0J  g  -  i  Jin  6  cosi?)"*  ***  d? 


Then 


where 


S  (vt'/u  «Cf> 

•im 


(5.4) 


(5.5a) 


(5.5b) 


«(Vl  = 


t aiQ+i.  8co»^ 


(5.5c) 


The  traveling  wave  nature  of  the  functions  can  be  seen  froa  their 

asymptotic  foras 


(5.6a) 


where 


I  2  r  z  -i* 

<=  V*V<«J  ftvU  ~  fartvfTi>w«f  J 


(5.6a)  results  froa  a  saddle  point  evaluation  of  (5.5a).  The  ratio 


♦  .  *L-  w  **»  « 


(5.6b) 


(5.7a) 


thus  converges  along  the  upper  branch  of  the  imaginary  axis,  da  v  >  0,  while 
the  ratio 


r-»  — 


:iv«. 


(5.7b) 


converges  along  the  lower  branch  of  the  imaginary  axis,  da  v  <  0.  The  product 
f.*Vi**»  hcwever ,  is  asyaptotically  not  of  exponential  character,  since 


r*  - ; — T 

b  f*  ~  mini 


(5.7c) 
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For  the  ratio  f'/-f£  in  (5.2),  we  can  then  obtain  the  series  representation 


V  $0 


Thus  different  expansions  are  required  in  the  upper  and  lower  half-planes. 
Multiplying  numerator  and  denominator  by  ,  we  obtain 


f' 

75  *  - 77775? - **> 


itrtV  £o 


It  will  be  shown  later  chat  a  stationary  point  of  the  Integral  J  in 
(5.2)  occurs  only  for  k  ■  0.  Thus ,  retaining  only  the  k  *■  0  tern  of  (9.8), 
and  introducing  the  integral  representations  (5.3a ,b)  and  (5.5b)  into  (5.2b), 
we  obtain  for  the  integral  & 

*4  =  T AWd-ijf  (5.9) 

where  }tj  represent  the  integrals  with  the  exponentials  respectively,  and 
....  ilUtVlP  % 

AW  *  ^ | 


*  it  *  $^-*<0**' +  **(n, ♦  «%-»)  J 

4.  =  '  (Vf’A)  1  V  (5.10) 

f,  *  i».  (ce*»  -i  *■>«  cats>x) 

$x* -(v,win(cos4  »i  smf  cotfx) 

* 

The  saddle  point  of  each  integral  is  determined  by  simultaneously 
equating  to  zero  che  partial  derivatives  of  the  exponent  S^2  with  respect  to 
each  of  che  integration  variables. 


_  as,  2  ^  _  sSe* 

iv  Sir,  i**x  if  Hi 


(5.11) 


to  determine  che  stationary  value  of  each  variable.  As  shown  by  Van  der  Pol 
and  Bresner  [Ref.  8]  ,  the  stationary  value  of  the  integral  is  given  by 


ii-jnr  c* 
h. I  =  AKx>  -755—  «  M 


(5.12) 


where  n  is  the  dimensionality  of  che  integral,  is  Che  value  of  che 
exponent  and  A^j)  the  amplitude  coefficient  at  Che  stationary  point  Vy , 
and  is  the  Hessian  determinant  [Ref.  10]  of  order  a,  evaluated  at  the 
corresponding  stationary  point 


*'*  *  |  ilX,  [  (5' 

From  (5.11)  we  obtain  for  che  scationary  point  of  Sj  ,  with  «fc»  v  >  0, 


(5.13) 


g  r  ;;v, ♦¥,)«. 
*’  =  X  (V,  *  *■*)*. 


(5.W) 
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and,  Introducing  the  stationary  values  ©i  *  *  0, 

w«%**i»*  8.»9  **  (5.15a) 

where  w1  and  v?  are  represented  as  w11 ,  wai ,  respectively,  to  denote 
association  with  the  stationary  point  of  Sj .  Froa  (5.15a),  wix  ,  wal ,  and 
(90+8)  forts  a  triangle.  This  leads  to  the  geooetrical  relation  shown  in 
Fig.  2(a).  Sinilarly,  for  the  stationary  point  of  S?  ,  with  ja  v  >  0,  we  oocain 
the  itae  relations  as  for  ^  ,  except  that  the  sign  of  8  is  changed,  so  that, 
denoting  the  corresponding  values  of  Wj  and  Wj  by  wia  and  waa  ,  respectively, 
we  have 


wit  ♦  "u  ♦  #»•  I  *  *  (5 . 15b) 

wla,  v33  ,  and  (8g- 8)  new  forts  a  triangle,  as  shown  in  Fig.  2(b). 

We  still  have  to  consider  the  portion  of  the  v-integration  for  which 
«fe  v  <  C.  For  this  situation,  we  need  cerely  change  the  sign  of  8o ,  correspond¬ 
ing  to  the  interchange  of  and  in  (5.8).  Then  (5.15a)  becoces 

w/  ■*•4  -  (®»-  8)  s  « 


or 


(5.15c) 


Similarly,  (5.15b)  beccoes 


(•R-w.i)  ♦  ( «-  w,',)  ♦  (»,(•»)**  (5. 15d) 

Froa  (5.15c)  and  (5.15b)  it  is  evident  that  w/j  is  the  suppleaent  of  wla ,  and 
w^x  is  the  suppleaent  of  w,a .  Siailarly,  froa  (5.l5d)  and  (5.15a),  it  is 
evident  chat  wla  and  va3  are  the  suppleaents  of  w^t  and  ,  respectively. 

Thus  the  geooetrical  interpretation  is  the  saae  as  for  «ia  v  >  0,  the  stationary 
points  merely  being  interchanged.  Consequently  the  result  is  just  twice  the 
contribution  due  to  the  stationary  points  in  the  range  As  v  >  0. 


The  relation  (5.15a)  for  Che  stationary  point  of  ^  can  be  fulfilled  only 
if  8,+  S  £  n.  For  9^+8  >  u,  Sy  does  not  have  a  stationary  point. 


The  values  of  Sj  and  and  v  at  the  stationary  points  are  then 


S)  5  D,  r*  j(8.*8> 
sj  =  D,t  ij(8.-e) 


I 

*  (5.16) 

j 


v,  =  Apy  vnw*  =  ipyvi«[-r- (S.  »«)-*.] 

if *^F7  U- iv »)-"«] 


) 


(5.17) 


where 


0,  *  -ttrtw,  -yccsv,. 


0j  *  -X  GGl  Wr,  -JtMWrt 
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Substituting  the  values  of  x  and  y  given  by  (2.21)  and  using  (2.7c),  we 
obtain 


0.*  -xktRce*  ♦  R.  cct  ;  -xk r, 

Oj*  -xk  (  3  cejw.j.  +  l?4  coi  wll)  t  -i.kr. 

For  the  far  field  we  allow  x  —  ®,  so  that,  as  can 


1 

be  seen  frocs 


(5.18) 
Fig.  2(b), 


r.  —  ft  *  R.tc*  tv-fl.-d)  =  R-RatasW.*®) 
r  _  R  ♦R4c*i(TT-6e»«>  =  R-R.ces(*a-6) 

»-(£,♦«)  w»»  -#*-  ‘6.-6) 

so  that 

v,  — •  (S.rt) 

v.  — *  -kR«v«(l.-3) 

Then  (5.18)  becomes 

o,  =  -xkfR-^cosiSsf  ei] 

Dj=  -*k[R_R.C6i(».-4)] 

the  only  non-zero  second  derivatives  of  S  required  to  evaluate  are 


1 

t 


1 


(5.17a) 


(5.19) 


the  follcwing: 

*>S. 

4*  Si 

- -  *  »CC1% 

TZf  =  »**•'» 

2TSU 

- =—  *  A 

*#• 

^=?i(«,tS«re 

*'S, 
— r  = 


Then,  denoting  the  determinant  &  by  |a.  | ,  where  the  order  of  the  elements 
i  and  j  is  v,  ,  u3  ,  q  ,  q, ,  v ,  ve  obtain 


VS*.  V.-.5 
Atx  vei.  ii«4 


1 


(5-20) 


Inserting  (5.16),  (5.17a)^  (5.20),  and  (5.18a)  into  (5.12),  using  the 
asymptotic  value  (5.7c)  for  f.ff*  la  A(vs),  including  the  factor  2  to  account 
for  the  contribution  of  Jm  v  <  0,  as  discussed  earlier,  and  introducing  the 
appropriate  value  of  (,  for  a  ring  source  (2rr  tines  the  value  given  in  (2.19)) 
we  obtain  for  E~  in  (5.1) 


(5.21) 


where  H(x)  is  the  Heaviside  step  function 


H  {<) 


*>  o 
*<0 


The  (unnormalized)  far-field  pattern  (i.e.,  the  dependence  of  Eg  on  3), 
thus  is  given  by 

'($}  *  (5-22) 

(5.21)  and  (5.22)  do  not  hold  in  the  vicinity  of  S  *  0,  since  then 
=  0,  corresponding  to  the  vanishing  of  the  second  derivatives  cPS^/dcl- 
In  this  region  a  higher  order  expansion  of  Sjp  is  required.  It  suffices  to 
evaluate  the  pattern  function  directly  for  3*0.  Then 

f  f**02  =  ?.-nW  =  • 

=  O 

Thus  there  is  a  null  along  the  cone  axis  for  the  ring  soiree  in  free  space. 

The  far-field  pattern  given  in  (5.22)  is  for  a  slot  of  inf initesimal 
radial  extent.  For  a  slot  of  finite  radial  extent,  it  is  necessary  to  inte¬ 
grate  the  field  across  the  slot  in  accordance  with  the  applied  field  distri¬ 
bution.  This  is  an  elementary  integration.  For  a  half-wave  slot  with  a 
sinusoidal  distribution,  the  result  is 

F(S)  *  Iwer*1  {tec  cos  [  |ccs(S.-«] 

(5-23) 

-a  to* l €,*« cc»£  j K (x-Se-Sj) 

The  first  term  becomes  indeterminate  when  3fl-  5  =  0,  the  second  when  5~"  3  “  tr. 
In  such  a  case,  the  indeterminate  coefficient  of  the  exponential  is  zero.  In 
the  neighborhood  of  x  =  0  (rr) ,  where  x  denotes  ?3-  5  (3d-  Smr)  .  the  critical 
coefficient  is  approximately  equal  to  's*/4. 


It  has  yet  to  be  shown  that  the  terms  in  the  expansion  (3.8)  for  k  >  0 
do  not  yield  stationary  points  of  the  integrals  Obviously  it  is  sufficient 

to  show  this  for  Ja  v  >  0.  3y  writing 


(<r*)v  = 


I 


using  (5.5b)  for  each  of  the  f{*'terms  in  the  numerator,  and  (5.7c)  for  the 
denominator,  we  obcain  for  the  exponentials  Sj^  in  (5.10) 


5ss  =  f  ter  *  £ 
»•* 

where 


4  s-iVf^JbolcsS^-iVrS.cos^) 
Then  cSj^/33,  “  0  leads  to  j,  *  0,  so  that 


41 


(5.24) 


(4J,  =  in+WBo 

Then  dS-  2/Sv  =  0  results  in 

wt+  (2<m)eo  i  a  =  if 

Obviously,  since  90  >  j,  (5. 24)  cannot  be  satisfied  for  k  >  0. 

The  evaluation  of  the  determinant  A,  given  by  (5.13),  can  be  simplified 
when  one  (or  more)  of  the  variables  is  uncoupled  to  the  others,  so  that,  for 
that  variable,  only  one  entry  appears  in  its  rcxv  and  column.  In  the  free* 
space  example  considered  above,  this  happens  for  the  variables  Oj  and  % . 

The  result  then  is  equivalent  to  inserting  the  asymptotic  expansions  for  the 
angular  functions  at  the  outset. 


5.3  RING  SOURCE  IN  SHEATH 

We  consider  next  a  ring  slot  on  a  cone  covered  by  a  sheath.  For  simplicity, 
we  consider  first  a  single  homogeneous  layer,  and  the  thin  sheath  approximation. 

For  the  far  field,  we  are  concerned  with  the  field  in  the  ambient  medium. 
Then  we  obtain  from  (2.11),  (3.7),  and  (3.15) 

"RE*  =  SO.iffA., lp<0  Arflpylvdv  I  S  ^xf 
«. /  ** 

*  5* <p<*} ^7)  ji„(p*)vdv  (5.25) 

®(v)  is  given  in  general  form  by  the  representation  (3.30).  The  first  term  of 
(3,30),  a,(v),  is  simply  the  free-space  term  (5.1)  (with  arguments  ox  and  py 
replacing  x  and  y,  respectively),  for  which  the  far  field  was  evaluated  in 
Sec.  5.2.  The  first  iteration,  Oi  (v),  is  given  by  (3.45).  On  inserting  this 
into  (5.25),  we  have 

-RE^  =(RRo)'/*M<  (5.26a) 


where 


d,  1 


IvlfxlK^py)  vdv 


(5.26b) 


The  integral  in  (5.26;*)  is  the  same  as  (5.2b),  except  for  the  replacement  of 
x  and  y  oy  px  and  py,  respectively.  Consequently,  by  comparison  of  (5.26a) 
with  (5.2a),  we  ran  immediately  write  for  the  field  El,  due  to  the  first 


iteration 


(5.27) 


where  E^,  is  the  frec-spaco  field.  Since  Roi?!  is  the  thickness,  t,  of  the  sheath 
at  the  source  radius, 


» 


(5.28) 


(j.27)  may  be  written  a3 
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(5.29) 


-y  (r.tfiE 


o 

1 


Thus  we  obtain 

•e,=  e;*e;,ou‘)- +  o(s*)  (5.30) 

where  Eq,  is  given  by  (5.21).  Thus,  in  first  approximation,  the  sheath 
decreases  the  far  field,  by  an  amount  which  is  proportional  to  6  and  to  the 
1  ’  square  of  the  electrical  thickness  of  the  sheath  at  the  source  position  on 
the  cone. 


5.4  SLOT  SOURCE  IN  SHEATH 


A  slot  source  covered  by  a  thin  sheath  was  treated  in  Sec.  3.3.2,  where 
(3.72)  and  (3.74)  for  c?(^)  and  cf (n)  ,  respectively,  were  obtained.  Whereas 
for  the  ring  source,  the  formulation  (2.11)  for  RI1D  involves  only  the  term 
n  =  0  because  of  the  syr.rsetrical  excitation,  for  the  slot  source  all  values 
of  n  are  involved.  Then,  as  pointed  out  following  (2.15)  and  (2.16),  the 
coefficients  qj3  and  qf  are  the  same  for  each  value  of  n.  The  relative 
excitations  of  the  various  orders,  n,  will  then  be  determined  by  the  applied 
excitation  as  a  function  of  the  azimuthal  coordinate  cp. 

Since  (3.72)  for  cf(u)  is  identical  with  (3.45)  for  the  ring  source, 
the  azimuthal  electric  far  field  of  magnetic  type,  E{2,  for  n  *  0  is  the  same 
as  for  the  ring  source  case,  which  is  given  by  (5.30).  For  n  >  0,  it  can  be 
seen  from  the  nch  order  asymptotic  expansion  of  (f'/fe), 

P£&ttot»)  ^  IjxntS11  *»«•-**♦ -nit/t) 

p^'(cos 8,)  U&»#  } 

^  jV  v< *•'*’  V  (->%  z  (5 . 3  la) 

that  it  is  only  necessary  to  multiply  the  coefficient  of  the  second  tern  in 
brackets  in  (5.21)  by  (~)n  and  the  entire  expression  by  cosSn*. 


For  the  slot  source,  there  is  also  a  meridional  component  of  magnetic 
type,  Ej?,  as  well  as  components  of  electric  type,  E^  and  E§.  Far-field 
expressions  for  these  components  will  now  be  obtained. 


For  Eg,  we  have 


ET  = 


i 


iW 


Rjrr#  36 


£.  r.  vnflC  f  . 


(5.32) 


The  Integral  differs  from  that  in  (5.25)  only  by  the  replacement  of  ^>'  by  jq  . 
The  asymptotic  expansion  of  is 

P„.*,(COS0'  /^inSA71  tO*(v6-7/*-iVT/z) 

P~"  ^(COSS,}  ~  * J,n* 

~  i[ e*^*.-*'  -  (-r; (5 .31b) 


which  differs  from  (5.31a)  principally  by  the  factor  v-1 .  Consequently  the 
saddle-point  evaluation  of  (5.32)  yields  for  the  nc^  term  (with  kg  =  k) 


t 


i 
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(5.33) 


The  field  components  of  electric  type  involve  o“(v).  q|(v)  is  zero,  so 

that  cpf(v),  which  is  0(6),  is  the  tern  of  lowest  order.  For  this  we  have  the 
integral  expression  (3.74).  Since  the  thin  sheath  approximation  is  limited 
to  terms  which  are  0(6),  we  need  be  concerned  with  the  evaluation  of  (3.74) 
only  to  terms  which  are  0(6).  Then,  as  shewn  in  Appendix  B,  (3.74)  reduces  to 


<nr  * 

For  Eg  we  then  have  from  (2.16),  (3.7),  and  (3.58b) 

c*  -  *  ftrc*  _  1  2nS,4. 6.  .«■ 


where 


in  which 


=  Z  i- 


=-i£^.  vh4Sd 


■  -»  7  h»Q) 


a,  =  2m  -  I  -  v 


dj=  2>n  1 1  -  v 
CLj=  -ZmM-V 


(5.34a) 


(5.34b) 

(5.34c) 


(5.34d) 


CU=  -Zm-l-V 

,  _  2 +v» 

n‘  "  »  *  “  Vlvf/*) 

*-n 

a  *  v(v-vi) 

h)*h4*0,  m  =  0 

h,  =  h2v 

/  m>0 

h<  =  h,/ 

(5.34c)  has  been  couched  in  such  a  form  that  the  poles  of  occur  in  the 

upper  half-plane,  Jn  A  >  0.  Furthermore,  since  it  turns  out  that  the  saddle 
points  of  (5.26)  in  the  A-plane  occur  at  the  poles,  the  lower  half  of  contour 
Cc  can  be  neglected. 

The  poles  of  W3V»±A  in  (5.34c)  pose  a  complication.  This  can  be  overcome 
in  the  following  way: 
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Denoting  the 


/. 

X-integral  by  I(80),  we  have 


HU. 


where  the  represent  the  poles  at  X  =  ±2m±l-v.  Multiplying  by  e*ak®o  and 
then  differentiating  with  respect  to  60,  we  obtain  the  differential  equation 


The  solution  of  this  differential  equation  is 

=  e'ia«8*  J  e*v  l.tyd^  *  f  ft ‘■°*t 


(5.35) 


where  the  lower  Unit  L  is  chosen  so  that  the  constant  of  integration  is  zero. 
Since  <Sm  =  -v  <  0,  ia^  has  a  positive  real  part.  Hence  x  ^wst  be  negative, 
and  L  =  -®  insures  that  the  integral  in  (5.30)  vanishes  at  the  lower  limit. 


Thus  (5.34c)  may  be  written  as 
*  =  y  t 

-rr  L, 

*•< 


where 


C  -  /]  «-**lu?4e^d<Te^Ve^  (5 . 36) 


In  this  form,  the  poles  have  been  eliminated  at  the  expense  of  an  additional 
exponential  integral,  which  increases  by  one  the  dimensionality  of  the  multi¬ 
dimensional  integral  to  be  evaluated. 


There  are  four  varieties  of  multi -dimensional  integrals  in  (5.36), 
corresponding  to  the  four  combinations  of  the  exponents  ±iv9  and  iiXSj .  Thus 
we  have  to  evaluate  the  four  integrals 

J1*.1*-  =  ftlm«di  (5.37a) 

1  ■€  I  • 

where 

ir.^-x  Zn'ii  m,#4= -z* *. I  (5.37b) 

Replacing  I, (px)  and  K^(py)  by  their  integral  representations  (5.3a,b)  as  in 
Sec.  5.2,  we  obtain  a  five -dimensional  exponential  integral.  By  equating  to 
zero  the  first  partial  derivatives  of  the  exponent  S  with  respect  to  the 
integration  variables,  we  obtain  the  following  equations  determining  the 
stationary  point: 


lift  r.’,  -  v/l(i 

(5.38a) 

vn  w.  =  a/»?y 

(5.38b) 

w.  ♦  *  « 

(5.3-Sc) 

=0 

(5 . 38d) 

v  =  >*  mk 

(5.38e) 
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On  inserting  (5.40a)  and  (5.40b)  into  (5.12)  to  evaluate  the  integrals 
we  obtain  for  E^  in  (5.34a) 


r«  -  ii.i  i 

~e  2t>  tir  4.  k'Vii  Jr 


1 ;  - . 1  £  ,» 


iin  Z.  —  * 


,h'  A-*.-**-'.*.-*-  _v_„ 


./  _  ,  ,  _  Os  •*!/  «  V* 

h. ,  -  v^rri  •  ■ '  _  : — i — : - r“ 

-1 


r*~o,  "-0-, 

In  a  similar  way,  we  obtain  for  e| 


i-;rs-  ; 


^>0 


E. 1 t. 


[•■  M- 

*X 


Wv'*  * 


<\-C^s.-9.-s)-  £;:'«,‘“i(2S*-s-'8)  *  ev-w'-e>1 


[v<-(2S4-e,«d>- H(T-e.-ri| 


(5.41) 


(5.42) 


The  far-field  expressions  (5.33),  (5.41),  and  (5.42),  as  well  as  the 
expression  corresponding  to  (5.21)  for  the  slot  source,  hold  for  a  slot  of 
infinitesimal  width,  corresponding  to  the  6(c)  factor  on  the  right-hand  side 
of  (2.13).  For  a  slot  of  finite  angular  width  23  <  X/2,  the  excitation  can 
be  assumed  to  be  unifom  in  c  across  the  slot.  Then  6(c)  in  (2.13)  is  first 
replaced  by  6(0-03),  where  03  represents  the  azimuthal  location  of  the 
infinitesimal  slot,  followed  by  an  integration  with  respect  to  c~,  over  the 
azimuthal  extent,  23,  of  the  slot.  Thus  i*  in  (2.19)  can  be  replaced  by 


The  factor 


2 


r. 


(5.43) 


thus  should  be  affixed  to  all  of  the  far-field  expressions  to  account  for  the 
O-dlstribution  of  a  slot  of  finite  width. 


5.5  EXTENSIONS  OF  THE  METHOD 

The  examples  worked  out  in  Sec.  5.3  and  Sec.  5.4  were  for  a  single -layered 
sheath,  and  in  Sec.  5.4  use  was  made  of  the  thin  sheath  approximation  to  0(6). 
There  is  no  inherent  difficulty,  however,  in  extending  the  treatment.  This  can 
be  done  to  any  order  in  6  without  invoking  the  thin  sheath  approximation,  since 
asymptotic  expansions  can  be  used  for  the  angular  functions  in  the  kernel 
functions  ^(v.X)  and  }^(v,X)  ot  the  single-layered  sheath  analysis,  or 
and  1 of  the  double -layered  sheath  analysis.  Thus  the  far  fields  and 
patterns  can  be  determined  for  the  general  case  to  any  desired  degree  of 
precision  0(6^). 
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SECTION  VI 


INPUT  ADMITTANCE  AND  MUTUAL  COUPLING;  EXTENSIONS  OF  THE  ANALYSIS 
6.1  INPUT  ADMITTANCE 

The  reaction  of  the  sheath  on  the  Input  adaittance  of  the  antenna  is  of 
great  inportance,  since  experinentally  It  Is  found  that  the  change  In  adait- 
tance  can  be  severe,  and  can  profoundly  affect  the  excitation  of  the  antenna. 
The  cethod  of  calculation  of  input  adaittance  was  presented  in  Reference  5, 
where,  for  a  radial  slot  energized  by  a  voltage  V  across  its  center,  the 
input  adaittance,  Y,  was  given  as 

Y=  -V*  Jdc  [  E,(6.)  6.dRo  (6.1) 

Since  %(So)  is  zero  for  the  electric-type  field,  this  reduces  to 

Y  »  -v*  Yd*  5  e;  16.)  h;«U  R.  si  nd.dk,,  (6 . 2) 

Since 

=  ZS±T,  *)  Rirxfj  (6 . 3) 

(6.2)  becoaes 

Y=  UwyV*)-*  (6-4> 

R!In  is  given  in  general  fora  by  (2.11),  and  for  the  input  adaittance  i  *  1. 

In  the  case  of  a  single -layered  sheath,  for  exaaple,  the  notation  of 
Section  III  is  applicable.  Then,  just  as  was  done  for  the  coefficient  Aj ,  the 
boundary  equations  oay  be  solved  for  Aj .  The  result  for  the  ring  source,  for 
exaaple,  can  be  expressed  in  the  fora 


A.M  = 

(v)  *  ^,(V)  +- 

(6.5) 

- 

(6.6a) 

= 

x  'Ss\ 

v  C,  if 

(6.6b) 

t-f)  = 

(6.6c) 

LpAc.A,)  : 

(6.6d) 

-0 

(6.6e) 

In  virtue  of  (6.6e),  tj (v)  is  0(6)  and  i. (v)  is  0(6n). 
Froo  (3.1a),  (6.5)  and  (6.6a),  we  have 


A 
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so  thac 


M*  =  [y.W ♦*»(»)♦  (6.7) 


Then  (6.2)  becomes 


Y=  Uu^T4  £.t.cci*s  E^U) R«  t.'. s. i R.- 
U''1-Y+){'i(.{V)  [*■)  t 


%+X«-Y*+ 


vhere 


X*5  (*"| ‘V,)'<  £,  £  t.cosfifdf5EIU,lR«*',‘i  **<*?,' 

<•8  •  » 

S(vx-W{^y(0^^^v 


(6.9a) 


(6.9b) 


*0  is  the  free-space  admittance,  so  thac  the  resaining  terns  represent  the 
effect  of  the  sheath.  In  (6.9a,b)  (f#/f;),  denotes 

/f.\  _  PJ^ulcostA 

P,^,»(cd‘6«) 

and  similarly  for  the  ocher  angular  functions  in  L(\  ,)a)  in  *-(v),  so  chat 
the  E  also  involves  che  Legendre  functions. 

There  is  no  difficulty  in  extending  the  analysis  to  the  case  of  a  slot 
source,  and  to  the  general  case  of  an  M- layered  sheath,  so  that  expressions  can 
be  obtained  for  che  input  admittance  for  these  cases  as  veil. 

Ec(90)  is  the  applied  field  at  the  slot  source.  For  a  single  half -wave 
slot  of  width  2v,  for  example,  ve  have  with  -}  *  v/R<. ,  where  Rc  is  the  locacion 
of  che  center  of  the  slot 


where 


E*  ^e*1=  '  Zw*  CcS k>(R.-^c'[ K{ ?«-RJ (6 . 10) 


Ri,t  *  Rc  *  7^/* 

w/Rc 


and  in  (2.19),  E,  3  -V/2v,  so  that 


f.=  -V/Hxik.Row) 


(6.11) 
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Since  t0 (v)  is  known,  all  of  the  integrals  contain  only  known  functions. 
They  can  be  evaluated  exactly  by  contour  integration  as  in  Section  III.  This 
procedure  is  not  very  appealing,  in  view  of  the  lengthy  series  representations. 
Since  the  admittance  is  strongly  affected  by  the  near  field,  asyaptotic  ex¬ 
pansions  of  the  angular  functions  are  not  suitable.  The  only  alternative 
appears  to  be  direct  numerical  evaluation  of  the  integrals.  In  this  connection, 
it  appears  chat  even  the  impedance  of  a  cone  slot  antenna  in  free  space  has  not 
been  evaluated  analytically. 

If  the  thin  sheath  approx imation  is  introduced  into  (6.6d)  for  L(>a,Xa), 
then  it  turns  out  chat  there  are  terns  which  are  0 (i^0),  and  these  cannot  be 
evaluated  in  closed  fora  either. 


6.2  MUTUAL  COUPLING 

The  coupling  between  a  transalccing  and  a  receiving  slot  on  the  cone  can 
be  expressed  in  terss  of  the  mutual  adaittance  between  the  two  antennas.  This 
was  discussed  in  Reference  S.  The  sutual  adaittance,  Yia ,  is  the  cocplex 
sutual  power  per  unit  voltages  across  the  two  slots 

Y,z=  Yn=  H  ( V,n  *  (V^  H  j)  dA  (6.12) 

Where  is  the  magnetic  field  strength  at  slot  1  produced  by  the  excitation 
of  slot  2.  For  example,  for  two  infinitesimal  slots  at  the  same  distance  Rg 
froa  the  cone  tip,  but  spaced  an  angle  i  around  the  cone 

(v,Vtr  J  B,,  (9.) H %&) RpiinBei 6,  (6-13) 


where 


E*  IBJ  =  E,  Kf)  +  Ex  i If  -  * )  (6 . 14) 

Thus  in  (2.13),  EL6(o)  is  to  be  replaced  by 


£,  Ho)  f  E*  S(f  -  §) 

Sy  putting  Ej  «  0,  it  can  be  seen  that  the  field  ^  produces  a  field  at  slot  1 
whose  ^-variation  is  shifted  by  i  relative  to  that  produced  by  a  field  applied 
to  slot  1.  Then  for  Hjj^Sg),  we  use  the  value  (6.6)  with  e?  replaced  by  o  -  J. 

Similarly,  for  two  infinitesimal  slots  on  Che  same  radial,  but  at 
distances  Rg  ,Rj  ,  respectively,  frca  the  cone  tip.  Eg  (55)6 (R -it,)  in  (2.13)  is 
to  be  replaced  by 


r  EzSfR-R.) 

Consequently,  with  Ej  “0,  it  can  be  seen  that  for  {^(Sg)  we  use  the  value 
(6.7),  in  which  the  i-  (v)  have  =  JtXr.Rj)  replaced  by  iSj/e.*) .*„(*, RJ. 

Thus  the  evaluation  of  the  mutual  admittance  between  antennas  can  be 
handled  in  a  manner  very  analogous  to  the  calculation  of  Che  self  admittance. 
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6.3  EXTENSIONS  0?  THE  ANALYSIS 

rniry  rhui- 

the  plasaa  has  been  idealized  bv  a ssuaing  that  it  <’<*  ?*Beral-  ln  tf>is  report, 
rad, a.  diction,  a*  cctinZT^lS^Ll viriac  a^L ”  “‘“'V*  *« 

steps,  thus  taking ^nt^afc^^^^ridio-^'^riati0  ffb£trar>‘  *****  of 
the  wave  equation  in  spherical  coordinate  ,ri  rlatio”*  ^  separability  of 
radial  direction.  The  radial  .arbitr4r>‘  variation  in  the 

variation  of  k  withi.  t£s  *  l  ^  £s  affected  ^  '»* 

tial  equation  has  to  be  replaced  bv  orj»°^bf 't/j3”  c°nstanc*  the  Bessel  differen- 
*  .  Then  the  K-L  traLfoL  aiso  cL  Z  tLZ  the  rad*a*  variation  of 

the  differential  equation  of  the  cylinder  *  <C  esploye,i‘  since  it  ste=s  fron 
applicable  to  the  new  radial  li  f  filial  ,  Unctions.  Instead,  a  transfora 
for  developing  such  a  transfora  exits  aL‘V*Uat£°n  15  required.  The  technique 
magnetic  Herra  v^ors !  “2”  “i“*  the  radial  electric  and 

this  difference  can  b«  ^lect^^*  STL  dlfferCnt£aI  equations, 

wavelength.  8  it  <  does  not  vary  appreciably  in  a  (local) 
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SECTION  VII 


CONCLUSIONS 

The  problem  of  radiation  from  slot  antennas  on  a  cone  in  the  presence  of 
an  inhomogeneous  sheath  has  been  solved.  In  this  repeat,  the  sheath  is 
considered  as  being  made  up  of  one  or  tvo  conical  layers,  each  of  which  is 
homogeneous.  In  the  method  of  analysis  used,  the  field  is  expressed  as  an 
integral  representation.  The  boundary  conditions  then  lead  to  a  system  of 
integral  equations,  which  number  4>f+A  for  a  sheath  composed  of  M(*  1  or  2) 
conical  layers.  By  an  extension  of  the  K-L  transform  technique,  these 
equations  are  reduced  to  singular  integral  equations  of  Cauchy  type.  An 
inversion  technique  is  developed  which  reduces  this  system  to  Fredholm 
equations,  which  can  be  solved  in  iterative  fashion.  3y  introducing  the 
parameter 

6.  =  1  -  o? 

where 

=  rt.x  /rt 

is  the  ratio  of  propagation  constants  of  adjacent  layers  of  the  sheath,  it  is 
shown  that  the  successive  iterations  proceed  in  powers  of  6,  .  For  a  suffi¬ 
ciently  fine  stratification  of  the  sheath,  the  first  iteration  should  suffice. 

In  general,  fields  of  both  magnetic  and  electric  types  are  generated  in 
the  presence  of  a  sheath,  even  though,  in  the  case  of  a  radial  slot,  only  a 
field  of  magnetic  type  is  generated  in  free  space.  The  field  of  electric  type 
then  is  created  at  the  sheath  boundary.  For  a  ring  slot,  hew  ever  in  which 
che  excitation  is  azicuthally  symmetrical,  only  a  field  of  magnetic  type  is 
,enerated  even  in  the  presence  of  a  sheath.  It  is  shown  that  the  solution  for 
this  case  forms  the  basis  of  the  solution  for  the  general  case. 

In  general,  the  evaluation  of  the  integrals  must  be  accomplished  by  contour 
integration,  which  leads  to  lengthy  series  expansions.  These  are  not  convenient 
for  numerical  evaluation.  For  the  case  of  thin  layers,  hew  ever ,  Taylor's  serlt- 
expassions  of  the  angular  functions  allow  ail  but  one  of  the  coefficients  to 
be  evaluated  in  closed  form. 

The  far  field  is  determined  by  a  suit i -dimensional  saddle  point  evaluation 
of  the  integral  representations.  This  is  illustrated  In  detail  for  the  free- 
s^>ace  case,  and  is  then  applied  to  determine  the  far  field  patterns  in  the 
presence  of  a  sheath.  This  can  be  carried  out  successfully  for  all  components, 
and  to  arbitrary  orders  of  iteration. 

The  calculation  of  input  admittance  and  mutual  coupling  between  trans¬ 
mitting  and  receiving  slots  on  the  con c  is  formulated  and  methods  of  calculation 
arc  discussed. 

Extensions  of  the  technique  to  more  general  situations  are  discussed. 
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APPENDIX  A 


THE  K-L  TRANSFORM 

The  K-L  transform  can  be  written  as 

.  A*  n 

f l*)=  —  A*  (AI) 

This  is  equivalent  to  the  transform  pair 

i U)  -  =  )  RvlIytxWdv  (A2a) 

F(v)  =  ?  (A2b) 

# 

Since  Ky(x)  is  an  even  function  of  v,  it  follows  from  (A2b)  that  F(v) 
likewise  is  an  even  function.  F(v)  must  be  an  analytic  function  of  v  in 
a  strip  of  finite  width  -6  s  Rev  s  6,  6  >  0. 

By  setting  f(x)  *  6(x-x0)  in  (A2b)  and  then  inserting  the  value  of  F(v) 
in  (A2a) ,  we  obtain 

.  *• 

i  (A3) 

A  second  6-function  relation  obtainable  from  (A2a,b)  is 

1  **(*>  **/*  *  l,,  WKvWdx/s  (A4) 

By  replacing  x  by  yR,  the  transform  pair  (A2)  is  frequently  written  in 
the  form 


In  terms 

f«9*^  J  Rv)1,(tR>v<v 

F(v)x  |  f(R)K„(V*}dMt 
of  the  spherical  functions 

(A5) 

(A5)  becomes 

xv(*R)«R* 

k,UR)*R'*  K,(fR) 

* 

i  *~T 

i(R)*5f  I  FMXvlrRJviy 
-*• 

(A6) 

and  (A3)  and  (A4)  become 

MR-R.)  =  JAiWR)-*,{rR,)vd» 

(A3a) 

-*m  -v*t 

(AAs) 

respectively. 


(A5)  is  equivalent 
[Ref.  6].  A  fact  which 
treatment  in  Ref.  6)  is 
the  proper  way  to  write 


to  the  form  originally  given  by  Kantorovich  and  Lebedav 
is  not  generally  realized  (again  steaming  from  the 
that  in  (AS)  y  is  to  be  considered  as  real.  Actually, 
the  second  equation  of  (A5)  is 


fm)  *  =  V‘pl'AAr,t\i*)K  v«*)4r 


(A5a) 


That  y  i*  be  considered  as  real  in  (AS)  is  forcibly  brought  out  by 
analyzing  in  detail  the  aanner  in  which  the  6-function  properties 

,0  R  SR* 


«... 

R» 

i  al  *,<«.<*» 

Ri 

are  displayed  by  the  integral  representation 

it 

R,  i(R-R.)  =  ^r-S  U*S)Ky(*R,)Vdv 


-1» 


(A7) 


(AS) 


In  (A8),  replace  Ky(Y8o)  by 


Then  the  integrand  has  the  asymptotic  behavior 

l.r,vl>o  (A9) 

The  firsc  term  of  (A9)  is  bounded  as  |v|  —  <*>  along  the  imaginary  axis.  If  je 
write  y  »  pei(P,  then  the  magnitude  of  the  second  term  of  (AS)  along  the 
imaginary  axis  is 


The  exponent  thus  is  positive  real  along  either  the  upper  or  lower  half  of  the 
imaginary  axis,  depending  on  the  sign  of  cp.  Thus,  the  second  term  of  (491  Is 
not  bounded  along  the  Imaginary  axis  unless  9  **  0;  i.e.,  unless  ^  is  real. 

If  <p  *  0,  the  magnitude  of  the  second  term  -*  Then  the  integrand  ol  (AS)  is 
bounded  along  the  imaginary  axis.  The  second  term  vanishes  on  an  infinite 
semicircle  .in  the  right  half-plane;  the  first  term  vanishes  aloog  an  infinite 

semicircle  in  the  right  half-plane  if  %  <  t0,  in  the  left  half-plane  if  R  >  Rq.. 

Since  the  Integrand  has  no  poles,  the  integral  vanishes  in  either  case. 

But  if  R  »  Rq  ,  the  first  term  of  the  Integrand  does  not  vanish  on  an 
infinite  semicircle  in  either  half-plane.  Consequently,  the  integral  is 
unbounded  for  R  ■  Rq .  But  if  an  integration  over  R  which  straddles  Rq  is  first 
performed,  a  factor  (v+l)-1  is  acquired  in  the  first  term  of  (A9)  which 

introduces  a  convergence  factor  as  well  as  a  pole  at  v  -  -1.  Then  the 

Integrand  vanishes  on  an  infinite  semicircle;  for  RS(>R0)  the  integral  in 
the  left  half -plane  encloses  the  pole,  and  it  is  easily  shown  that  the 
residue  is  Rq.  ,  The  seme  result  is  obtained  if  an  integration  is  performed 
over  Rq  over  a  range  which  straddles  R.  Thus  the  6 -function  properties  (A7) 
are  ail  properly  exhibited  if  y  is  real. 

The  integration  over  R  over  a  range  straddling  Rq  ,  as  sketched  above, 
also  provides  the  means  whereby  the  assumption  chat  y  is  r -.al  is  justified. 
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It  is  merely  necessary  to  integrate  R  along  a  contour  rotated  by  the 
angle  -©,  as  in  (A5a) ,  which  makes  yR  real.  This  is  merely  another  aspect  of 
the  6-function  property  of  (A3);  i.e.,  it  does  not  make  "sense”  mathematically 
until  an  integration  ever  R  (or  Rq)  is  performed.  In  the  case  of  the  K-L 
transform  of  a  general  f(R),  the  justification  for  considering  y  as  real  is 
performed  by  the  inverse  transform  (A5a) . 
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» 


SOLUTION  OF  THE  BASIC  INTEGRAL  EQUATION 

The  basic  integral  equation  is  of  the  fora 

J  A{V!  -«*i/>vdv  «  S  B(>;J  i* (px)ydv  (Bl'j 

where  A(v)  and  B(v)  are  even  functions  of  v,  and  contour  C-  is  the  imaginary 
axis  of  the  v -plane.  It  will  be  assumed  that  p  is  real,  in  order  to  avoid 
convergence  difficulties  in  later  inversions  of  the  integral  equation  derived 
fron  (Bl).  This  can  be  justified  in  the  case  way  as  in  Appendix  A,  by  noting 
that  eventually  integration  over  the  source  coordinate  will  be  required. 

Hence  in  the  ultioate  integration  over  the  source,  the  oath  need  merely  be 
rotated  so  that  px  is  real. 


iv(px)  in  the  right-hand  integral  of  (Bl)  is  now  expanded  [Ref.  7)  as  the 
series  in  iy*2«(x): 


=  **  w  =  x-VZ  c_(v,p)  CWMWUi 

where  the  coefficients  c,{v,p)  are  polynomials  in  pa 

C«(V,p)  =  K- 

<e!r(vtii  1,1  ‘  ‘r>  fz  ni+trji  (o-jj*  jt 


(B2) 

<E3) 


•pie  following  properties  of  these  coefficients  are  important  for  later 
deve lOpments : 


By  writing 


6  *  I  -  p* 

(B3)  may  be  expressed  in  the  form 

£ctv)Sr 


(B4) 


(B3a) 


where 


*  «•!  )Tv»i»»j 

c.(v>  =  cj(vl  *  « 
cl  tv)  r  O  \ 


(B3b) 


cj.lv)  «  -V 


«>0 


Thus  the  first  ter*  It,  the  expansion  (B2)  is  0(1),  while  the  remaining  ter 
are  0(6). 


With  (B2),  (Bl)  beccoes 


/ ********  -  if'  *v>  Murti  w  ix  (35) 

^  C*  <U|  ^ 

Taking  the  K-L  transform  of  (B5) ,  by  multiplying  by  V(x)d x/x?  and  integrating 
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over  0  to  «,  we  get 


AlpJs^  j,  (B6) 

wi.'jte  u  Is  i.  -jr  t  '*ai  ~o{.ot  on  i-  ■>  l*^jsir.*ry  axis ,  Cj  is  5  contour  parallel 
to  *.U  M'a'rit.  .ry  axis  tc  ‘he  rij  .  oC  #,  as  shown;  '  '  FU'.  Bl. 


Fig.  Bl 


and  H(v,p;p)  is  given  by 


X(u  i-t)  =  £  ?  Sm£*m£. 


(B7) 


Alternatively,  instead  of  expanding  iy(px)  in  (Bl)  in  terms  of  iy^s, (x) ,  we 
can  expand  i,(x)  in  terms  of  c  series  in  iv.»2a(px)»  This  can  be  done  by  replacing 
x  in  (Bl)  by  y/p  and  using  (B2) : 


Then  (Bl)  becomes 


=  *  fwE#c* —3-  iy„„ty) 


(B8) 


IBM i, (,J v SV  -  f  fwAM  £  C.tV,V|) (V4l»)  dv  (B9) 

•** 

The  transform  of  (B9)  then  is 


[AW  X(mv;  WdV  (BIO) 

Thu*  (BlOl  it  she.  Inversion  of  (B6).  and  vice  versa.  X(p ,v; 1/p)  is  given 
by  (B7)  with  the  replacement  of  p  by  1/p, 

It  should  be  noted  chat  X(p,v;p)  is  e<{uivalent  to  the  integral 

*  vXiy(fs)  JkfWWt1  (Bll) 


Similarly, 

ft)  =  vf  4»/*»  (B12) 

In  Sec.  3.3,  tue  related  function  is  encountered.  This  is  defined 

"  ?  5  ^C|iel,v;f)  ♦  ^  X<H-f:viPl]} 


by 


Using 
che  a 
a  *  0 


(B7)  for  the  X-funetions,  ic  is  obvious  froa  (B3a)  and  (B3b)  that  only 
*  0  terms  are  0(1),  the  remaining  terns  being  0(6).  Separating  off  the 
terns,  we  have,  since  c0(v,p)  *  i. 


J _ \ 

V-I 


I 


_  (_ l. 

?  I  H  * 


v  ♦,*  f 


.‘Us-snsffl 


(B13) 


In  a  typical  integral 


S  r(v)3t  dv 

where  F(v)  is  an  even  function,  we  can  change  the  sign  of  v  in  the  second 
tern  of  each  pair  of  terns  in  large  parentheses  in  (B13),  for  which  the  contour 
becoaes  C/-1  (see  Fig.  31).  A  shift  of  contour  froa  c/-l  to  Cj+1  allows  these 
terns  to  be  coabined  with  the  first  tern  of  each  pair  in  (B13) ,  and  also 
collects  residues  at  the  poles  v  «  Liu •  The  residues  cancel,  leaving 


S  FWftf***'! 
c*.* 


-r~>-  ? 


(B14) 


Since 


f*-pw=  ♦<>«*> 

each  of  the  terns  of  the  integrand  of  (B14)  is  0(6).  Hence  integrals  containing 
K±U.±-J  are  0(5). 

In  subsequent  developments,  we  encounter  double  integrals  of  the  type 

dtp)  *  —  S  THrrfk}*  (*>5) 

where  M(v,X)  is  an  even  function  of  both  v  and  X.  Since  it  has  just  bees 
shown  that  the  inner  integration  yields  a  result  which  is  0(6),  this  result  may 
be  represented  as 

IEiv> 


where  E(vy  is  an  even  function  of  w.  Hence,  if  we  wish  to  evaluate  **&)  to 
0(6) ,  ve  need  be  concerned  only  with  the  first  term  of  the  series  representing 
X(ji,v;i/p),  since  the  coefficients  of  all  higher  terms  of  the  series  are  ail 
0(6)  (see  (B3b)).  Thu*  we  have 


Since 


ve  find 


=  55X  r*«v»  *'  ♦  otw 


fvs  ,-!><  »0t« 


*w£[S-5)TT*i*o«*  -*e«p»+°c« 

^  f'r 
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Consequently,  to  0(6)  (B15)  cay  be  written  as 

=  t  0«*J 

Similarly,  for  integrals  of  the  type 

^  *  "ST 

S  ci 

where  ^(v.v)  *  0,  so  that  the  inner  integral  is  0(6),  we  obtain 

t  0(S *) 


APPENDIX  C 
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APrtKDIX  D 


ZEROS  0?  HE  LECEXDXE  FUNCTION* 


The  purpose  of  this  progran  is  to  confute  the  zeros,  v, ,  of  the  associated 

Legendre  function  of  the  first  kind.  PY~(c .  as  a  function  of  degree  v, 

and  tie  derivative  ^-?T~(cos3)  at  v, .  A  method  of  computing  these  *era» 
was  described  by  Vile ox*  of  the  University  of  Michigan.  However,  a  different 
method  is  employed  here,  which  is  acre  accurate  and  faster  when  9  is  close 

to  180  degrees,  which  is  the  region  of  interest  here.  Double  precision 

arithmetic  is  used  in  order  to  achieve  the  desired  accuracy. 


The  following  equation  was  used  tc  evaluate  the  associated  Legendre 
function: 

a  r( *•««-✓, r,  y.  ,, a,  .  , (oi) 

rt-rz  z — 

Then  os3)  is  given  by: 

6)  =  ^cct JTkcsi) - i  */-. vs  5-f  “***§ 


I  ’  ,"rt»-»)r{»»ln5«!(r.*i4!  ,l  *  *  *  I 


(02) 


In  these  equations, 

3,  s  *yfr-»)-  *(»«*l*v)  (03) 

>,  *  -  yfenei-vl  ♦f(o-V)-  (I>;) 

k«  =  ’>(«»;!  ♦  •f<r.t9t»0 -♦{«•*««- s)  -V,{«*s»sUv)  CD5) 

V.'  5  «■. -  Vl  -  v/(r.  +  <-.  ►<  fv)  (Do) 


The  first  tens  In  the  derivative  vanishes  at  v5  ,  of  course.  Bcwever,  the 
complete  expression  is  needed  in  the  iteration  process  to  refine  successive 
trial  values  of  the  zero. 


In  order  to  take  advantage  of  the  expressions  which  are  common  to 
?^*(ce*S)  sad  JLp^HcosS'  .  the  arrays  ZXA,  123 ,  ZXL  are  set  up.  These  are 
defined  as  follows: 


ncwjr{jv.uy? 


f 


ZR*. 


rfe»<»-yi 

r  irr^i  Tln+iri  fil  »«Jt 


U» 


U«>.S  i 


ZXL,  =  | 


In  terns  of  these  quantities. 


(D7) 

(08) 

(09) 


f' 


* 


*  FETES  8.  WILCOX,  "The  Zeros  of  F$  (cos8)  and  0/SS)£  (cos5) ." 
Mathematics  of  Computation,  Vol.  22,  So.  101,  Jan  196S. 


* 
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(5jm(cos6)  =  s>n vn  tafl'f  [L‘ ZRA.  ZRB.ZKU]  (DIO) 


and 


£  P;"’(C0»8)  =  -ir  V«vtr 

+  T  COt  ^•n  P„  ’'(cosfii 


tan  *  f  [l  ZRAp  g.»  (.)-  £ZRS„  (ZKL.  h.  f  o] 


(Dll) 


The  ratios  of  the  gamma  functions  in  ZRA  and  ZRB  are  computed  using  the 
recurrence  relationship,  so  that 


r(rt-v)  r(n«-l  »v)  /=?  ... 

r(m-v)ri»->.uv)  =  /JJ  (w-J-vKm-iti**) 

and 


(D12) 


ftn+r*-*)  rlnHrl»l  W) 
r(m-v)  riintiTJ) 


TT  (m+itii)  (m+i-l-%)) 


(D13) 


Subroutine  XP  computes  P^®(cos6).  and  in  the  process  computes  the  arrays 
ZRA,  ZRB,  and  ZKL  by  calling  other  subroutines.  XP  calls  three  other  subroutines, 
RA  to  compute  the  ZRA  array,  RB  to  compute  the  ZRB  array,  and  EK  to  compute  kn, 
which  is  needed  for  the  ZKL  array. 


Subroutine  XPPR  computes  A-Pjla(cos9)  using  the  arrays  ZRA,  ZRB,  and  ZKL. 
It  calls  three  other  subroutine^,  EKP  to  compute  k„,  G  to  compute  gn,  and 
H  to  compute  hn. 


The  EK,  G,  and  H  subroutines,  which  compute  kn,  gn»  and  hn,  respectively, 
all  call  the  PSA  subroutine,  which  computes  the  f  function.  The  EKP  function, 
which  computes  kn,  calls  the  PSB  subroutine  to  obtain  the  t'  function. 


The  if  function  is  computed  by  using  the  recursion  formula  to  step  tne 
argument  upward  until  it  is  greater  than  10,  and  then  using  the  asymptotic 
formula.  Thus  the  formula  used  is 


=Sfl3(x^)  -  £  -  ^ - £(**»)' 

xtn  >10,  =  &2m/2m 

where  the  Bg,  are  the  Bernoulli  numbers. 

$ 

Similarly,  the  if'  function  is  computed  using  the  formula 

V'M  =  zb  (tb5  t*"4  (*«♦!)»  +  illb?  B,;(z»«)',<-'' 

Z*n  >  10,  B*  =  Zw  Bm  =  Bjm 


(D1A) 


(D15) 


To  compute  the  derivatives  with  respect  to  0,  the  recurrence  relations 
are  used. 

P,""(cosJ>  =  (D16) 

and 
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where 


~  P^Uoit)  *  TfT)  [h  frpccD-tu-rhjO'"*! +  cM6?;"^i6)-f^Wi]  (D17) 
In  the  actual  computations,  the  first  values  computed  are 
jf  P^(Wi9)/(-A£^3)  aT^Pr(t«3)/(-A8/ii«6) 

-AB  =  -  (tan"  i  S'>V*)/r r 


Then  in  (D16)  the  following  is  used 


; lv-m)  P‘,  (o» 6)  - 


I 


St«(v-I)ir 


(v-ip) 


and  similarly  for  the  corresponding  term  in  (D17). 


The  program  consists  of  a  main  program  and  ten  subroutines.  All  routines 
are  written  in  Fortran  IV. 


The  main  program  first  reads  a  single  control  card  which  contains  the 
parameters  for  the  first  set  of  zeros.  There  will  be  one  additional  control 
card  for  each  additional  set  of  zeros.  The  program  terminates  by  reading  an 
end  of  file  while  attempting  to  read  another  control  card.  A  description  of 
this  control  card  is  given  on  page  67. 


The  zeros  of  Pj^cosS)  are  denoted  in  the  program  by  Nu,  and  the  zeros  of 
JLpJ^cosB)  by  Mu.  The  program  uses  an  initial  (estimated)  value  of  Nu  from 

the  control  card  to  compute  the  first  values  of  P^fcosB)  and  ^^(cosB). 

The  derivative  is  then  used  in  the  following  way  to  obtain  a  second  value  of 
Nu  which  gives  a  Py^cosQ)  with  a  smaller  absolute  value: 


=  <  ~  /(  i*  P»  h*  v,., 


(D18) 


This  process  is  continued  until  P17n(cos8)  is  sufficiently  small.  The  criterion 
of  smallness  is  selectable  by  means  of  a  tolerance,  which  has  been  set  at  Iff*10  . 
An  additional  tolerance  of  Iff*8  for  the  second  term  in  (D18)  is  also  Included. 
This  insures  that  8-decimal  accuracy  is  obtained  for  the  roots.  These  tolerances 
are  arbitrary  and  could  be  changed  if  a  different  accuracy  is  desired.  Usually 
only  a  few  iterations  are  required  to  obtain  ten  significant  figures  of  accuracy. 


In  a  similar  way,  using  an  initial  value  of  Mu  supplied  by  the  control 
card,  a  value  of  Mu  is  found  for  which  ^JU1J‘(C0S8)  is  sufficiently  small. 

The  four  values  for  the  fiist  zero,  v0 ,  ^-Pv“?,  an<*  are 

O V  qVOw 

stored,  and  the  next  zero  is  then  processed.  The  initial  values  of  Nu  and  Mu 
for  the  next  zero  are  obtained  from  the  final  values  of  Nu  and  Mu  for  the  zero 
which  has  just  been  processed,  tt/6  is  used  as  an  approximation  for  the 
difference  between  consecutive  zeros  for  both  Nu  and  Mu. 


When  all  the  zeros  have  been  processed  chat  were  called  for  in  the  control 
card,  a  page  of  output  is  created  containing  a  tabular  listing  of  the  four 
values  corresponding  to  each  order  of  zero.  The  program  then  attempts  to  read 
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a  new  control  card.  This  is  the  usual 
various  stages  of  the  computation  are 
entries  in  the  control  card. 

A  flow  chart  of  the  urogram  is  on 
listing  on  pages  69  -  73. 


procedure.  Additional  printouts  at 
selectable  by  means  of  appropriate 

page  68,  and  the  complete  program 


The  control  card  is  read  by  the  following  Fortran  stateaent: 


READ  (5,55)  KW,NZ,M,AA,BB,CC 
55  FORMAT  (5I1.2I5.3F10.5) 

KW  is  a  dimensioned  variable  (DIMENSION  KW(5)> 


The  parameters  in  the  control  card  are  described  below. 


Col.  1  Non-zero  meats  print  values  of  v,  P^®(cos0)/(-A3) ,  and 

—*  Fy  (cos8)/(-AB)  for  each  trial  value  of  v  for  each  zero, 
dv 

Col.  2  Non-zero  means  print  final  values  of  v,  P^"*(cos0)/(-AB) ,  and 

P^*(cos8)  for  each  zero. 

3v 

Col.  3  Non-zero  meant  print  values  of  p,  l£*(cos8)/(-AB/sin9) ,  and 

OO 

v2 

-----  ^7*(co*9)/(“AB/sin0)  for  each  trial  value  of  p  for  each  zero. 
opo8 

Col  4  Non-zero  means  print  final  values  of  p,  l£*(cos8)/(-AB/sin8),  and 

00 

JL.  aT*(cos0)  for  each  zero. 
ojxoB 

Col.  5  Non-zero  means  print  a  page  containing  final  values  of  v, 

P^^cosS),  p,  -3* .  ^*(cos8)  for  each  zero  in  the  set.  This  is 

dv  3p38 

the  normal  production  output. 

Col.  6-10  Number  of  zeros  in  the  set  to  be  computed. 

Col.  11-15  Order  of  Legendre  function  (m  in  P^^coeS)). 

Col.  16-25  Argument  in  degrees  (8  in  P^*(co*9)).  This  should  be  greater  Chan 
160  and  less  than  180. 

Col.  26-35  First  trial  value  for  v  in  first  zero  of  the  set. 


Col.  3 


Col  4 


Col.  5 


Col.  36-45  First  trial  value  for  p  in  first  zero  of  the  set. 
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t  JOB  178l01-A,25,lO(5O,JO*N$ 

BEX E CUTE  _  1BJG8  _  _ 

BIBJO'B '  ' 

BIBFTC  TEST 

ooubls  precision  Pi,ANC,THtENU,£HU.A,8,s.sr' 

00U8LE  PRECISION  Q.0PA.2RAf2RB,ZKL« XPNU.XPNU 
oOUBLf  PRECISION  C,CT,fNtCSO.TOL,P,PA,PFRtPFRA,W 
_  01  PENSION  2RA(10I,ZR_B(IOOI,ZKL(100) 

DIMENSION  KM (57 

OINENSIGW  XNUC 200 ) »XKUI200I ,XPNU(200 )»XPRU(200) 
777  CONTINUE 

RE AO  (9*53)  KM,N2,M, AA,8B,CC 
33  FORMAT  (5tl«213»3F10.3) 

EH-M  _ _ _ 

*  ANG-AA  "  '  "  . 

ENU>9B 

EMU-CC 

TOL-l.OE-lO 

TT0L-1.0E-8 

PI-2. 1*15926 33 389793  _  _  _  _ 

"  TH-ANG«Pl/28~G. 

CT-OCOS(TH) 

ST-OSIM(TH) 

C-OCOS(TH/2.) 

S-OSIMCTH/2. ) 

CSO«C*C„  _  _ 

'  00  300  J-l.NZ 
N«i 

1  CONTINUE 
A-OSINCENUMI) 

CALI  XP(NfENU«CS0»P»2RA»2Ri»2KLtNNI 
CAU  XPPRIN*  ENUt CSQ»PPR»2RA»ZRBr2KL»NN) 

PPR-PI  •P*OCOS  flNU*P  if/AiPHK 
IF  (Ktfil).EO.O!  GO  TO  11 
MRITE  (6,44)  N,ENU,P,PPR 

11  CONTINUE 

IP  (ABS(P).LT.TOL)  CO  TO  888 

IP  (ABS(P/PPR) «LT .TTOLI  GO_TO_  888  _  _ 

N-N«l 

IP  (N.CT.10)  CO  TO  888 
ENU-ENU-P/PPR 
GO  TO  1 
888  CONTINUE 

8-((S/CJ**MJ/PI _  _ _ 

»R--MPPiT 
PPR»PPR*A 
XNUiJI-ENU 
XPNUI JI-PPR 

IP  (KMUI.EC.O)  CO  TO  12 

UNITE  (4f (i4)_R *ENUfP»PjPR _  _ _ 

44  FORMAT  ( 110, 3C 30. 141 

12  CONTINUE 
N-l 

2  CONTINUE 
A-OSIN(ENUPPI) 

__  CALL  XP(N«Ej|UjCSQ.f Pf ZRAf ZRBtZKLyNN)  _ 

CALL  XPPR(N*EMU»CSQtPPRf 2RA»2R8«2KUNN) 

CALL  XP(M, EMU-1., CSQ,PA,ZRA,ZRB,ZKL,NN) 

CALL  XPPRfM, EMU-l.nCSGfPPRA, 2RA* 2RB«2KL*NN) 
Q-EMU*CT4P*(EMU-EM»*PA 


T 
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Q^-EMUKTMF*4(£MU-E>04FFA'A#Cf*FVFA 
QgMQ»f!j>pCp$  ( INU*FI  )/A«QFR  _  _  __  _ 

if  o^m.EoVo)  co‘  ta'iS' 

MUTE  *56,4*1  N,ENU,Q,OFR 

'  13  cmivtir 

IJ^(A6:*CQJ,IT.T3L)  CO  TO  499 
IF  (A».ftQ/OFM) •LT.TTOLI  i Jff  f<5  *49 

N«N4l _ _ 

IF  (N.CT.IO)  CO  TO  444  . . 

ENU-ENU-Q/QF* 

GO  TO  2 
949  CONTINUE 

OFR— MOFR/S? 

QFMQFFM _ _  _ _  __ 

XNUIJI-ENU 
_XFMUUI-QFA 
IF  (KM(4).EQ.O)  CO  TO  T* 

_ MITE  (4,441  N,ENU,Q,0FA 

14  CONTINUE  ' 

£WENU*H/TH _ 

enu-enuffi/th 

300  CONTINUE  _ _  _ _ _ 

IF  (KMISI.EQ.O)  CO  TO  IT  - 

MUTE  44,1411  44, N _ _ _ 

141  FORMAT  ( 1H1 , 2Q  X, 20HAN6LE  IN  DECREES  IS  ,F9.9»20X,«H0R0ER  IS  ,TsT 

MITE  (4,1421 _ _ 

142  FORMAT  (/9X,3HZNO,20X,2HNU,20X,4HF«NU,20X,2HNU,?9X,4HF*NU/71 
MITE  (4,1431  (I«XNU(  I  I,XFNU(  1 1  ,XNU(  I  )  ,  XFNUf 1 1 , IU,NZ » 

143  FORMAT  1 3X,1 3, 15X,FH.T, I IX, 014.?,  1  OX.F l 1 . ?,ll X . D 1 4 . 7 » 

MITE  (4, 1441  _  _ _  _ 

144  FORMAT  (INI)  . . 

MUTE(12.1T» I  NI,M,AA _ _ 

170  FORMAT  (2I9,F10.9)  * 

MITE(  12,177)  (!,XNU(!i,XFNUm,XMO(I),IFWUm,?-l,NZ> 

177  FORMAT  (I3,4017«10) 

19  CONTINUE 
60  TO  777 

_ PW _ 

•IOFTC  XF. 

SUM OUT I NE  XF(M,ENU,CSOfF,ZRA,ZRR,ZKL,NNi 

OOUfUE  FRiCtSlCN  F(CSQ,CNU,RA,lRA,Zt,RT7MI,ER,Za(C7r7i^l 

OOUIUE  FRICISION  TOt  _ 

01  Mi  NS  I  ON  ZRAliO)  ,ZRR(100)~,ZKL(100) 

TQt«l,Qg«10 _ 

AM. 

IF  (N.IT.1)  60  TO  101 
00  100  NN«i,M 
Urnm-l 

ZRA(MI)*RA(N,M,TNU,CS0) 

A»A,ZRA(NN) _ _ _ 

Iff  CONTINUE 

[01  CONTINUE  _ _ _ 

MO. 

NN«1 _ _  _ _  _ _ 

300  CONTINUE 
NbNM-1 

ZR8(MN)aRlfN«M«ENU,CSQ) 

ZKL(NN)»EK(N,F,ENU)-0i06(CSQ) _ 

Z0«ZRt(NN)4ZKL(NN) 

»>I*ZI 


I- 

i 


IF  (NN.GT.99)  GO  TQ  201 

IF  (A6S(Z8).LT.T0U  GO  TO  201 

NN-hVi 

_  _  GO  TO  300 

201  CONTINUE 

WRITE  (6, *6)  NN 

66  FORMAT  (20K  NUMBER  OF  TERNS  IS  * 14) 

_ 88«(-ll»»M  _  _ 

P*A»BB*B 

RETURN 

END 

»I8FTC_XFFR. 

SUBROUTINE  xFFR(NeENUfCS0;F>R,TRA,2RB,2KL#LU 

_ DOUBLE  FR£CISICNL_CSQ*ENU,_FFR,ZRA,2RB,2lCL.ZH,ZEICF»EI(F,_G,H,BtCA1C8 

DIMENSION  ZRA (10)  *ZR8(  100)  *ZKL(i30~) 

8*0. 

IF  (M.LT.l)  GC  TO  101 
00  100  NN-1,M 
N-NM-1 

8«8«ZRA(NN)»C(N*M*ENU> 

100  CONTINUE  "  ~ 

101  CONTINUE 
CA* (-1 )*»N 

cb-o. 

00  200  NN«1*LL 

N»MN-1 _ _ 

ZEKP*EKP(N»N,ENUJ 

ZH*H(N*M* EMU) 

'  CB-CB»ZRBt<fM>«(ZKL(NNf*ZH*ZEKFf 
_200  CONTINUE 

'ppr-bVCa*Cb 

RETURN 

!(«'•  ’  "  ' 

SIJFTCRA. 

DOUBLE  PRECISION  FUNCTION  RA(N*N,£NU.CSO I 

_ DOUBLE  PRECISION  RAy£NU*CSQ*PROO*EJ»El*EMvEL 

DOUBLE  PRECISION  FACT fill 

_ DATA  (FACTII )*I»1*11)/1. 00*1.00*2.00*8.00*29.00*  __ 

1  12.01 * 72. 0 1*504^01*40  32 .0 1 *  38288.01* 38288.02/ 

EN-N 

K«M“N  —  . .  "" 

PROC-1. 

IF  (K.lTTI  )  GO  TO TOi' 

00  100  L«1*K _ 

EL-L 

_ PROO«FROO/f  tfF^EL-fNUI8lEN-ELU.J>EMUI  1 

100  CONTINUE 

101C0NTINUE  _  _ 

F*CD-FROO*WCT  (RJ/FaCTTmT) 

_ FRCO"FR009U"CSOI»*N) _ 

RAMROD  * 

RETURN  _ _ _ _ 

ENO 

8I8FTC  R». _ _ _  __ 

OOUILE  FRECISICN  FUNCTION”  RBlNfNiHlU.tSfll 

DOUBLE  PRECISION  RB.EWU.CSQ.A.EK. EJ.EL, EN _ 

OQUBLE  PRECISION  FACTIlil 

DATA  fPACT(I)*!-l*U)/l.OO*1.00j>2.00»8.00*24.00* 

1  12.01*72.01 *504.01 *40 32. 01 *38288.01* 38288.02/ 


IF  (N.6T.0I  00  TO  1 


RB-CSQ**M/fACTtX*l) 

GO  TO  999 

1  EN-N . .  -  - 

__ej-**n 

EL«F*N-1 
_A-EN»EJ 

RS-RB*CS0*TE  j*  CNU  )•  i  Et-EHU  J/A 

999  RETURN  _  _ 

ENO 

SIOFTCJSA. 

OQUBLE  PRECISION  FUNCTION  PSACZI 
OOU8LE  PRECISION  2*X,B,SUN*PSA 
'  '"DIMENSION  B€  31 

DATA  (BU), 1-1,5)/ 
l  0.B33333A333 333333  0-0 l, 

2_-0. 8333333333 333333  0-02* 

'3  0.39*82539*12539*1  C-02, 

_ * -0.41****&**t*****7  0-02* 

5  0.757575757575757*  0-02/ 

X«1 _ _ _ _ 

SUMO.  • 

_1_IF_CJU6T.10._I  CO  TO  10 

sum-sun-Y./x  -  ~ 

X-X*l. 

CO  TO  1 

10  SUM-SUM*DL06(X)-l./tX«X} 

DO  100  1-1,3 
11*1*1 

“  “SUM-sun-»<i»/Tx**TiT 

100  CONTINUE  _ 

PSA-C-UN  . . 

R27URM _ 

(JO) 

OOUBLE  PRECI  STON  PuNCTTOnPSSU) 
DOUBLE  PRECISION  I, X,8,SUH,PSB 
DIMENSION  8(5 1 

_ OATA  («f  11*1-1.51/ _ _ 

1  1. *************** 7  0-01, 

2  -0.3333313311333333  0-01,  _  _ 

3  0 . 23809 5 23 80952311  0-01, 

*  -0.3333333333333333  0-01 , 

3  0 .737575737573757*  0-01/ . 

X— Z _ 

SUM-0. 

__  L  IE  f  X.CT.1Q.  1  60  TO  JO 
SUM-SUM  l./TftX) 

X»X*1. 

GO  T31  '  '  "  " 

10  SUM«SUM*1./X*.3/(X*XI _ 

00  100  1-1,3 

_ II- 1*1*1  _ 

SUM- SUIUB  0171 X**  Iff  ' 

100  CONTINUE  _ 

PSB-SUM 

_ RETURN _ 

ENO 

81BFTC  IK* _  _ _  _  _ 

DOUBLE  PRECISION' FUNCTION  EX(N,M,ENUT 
DOUBLE  PRECISION  PSA 
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OOUBLE  PRECISION  ENU,E*»EN,EH,ZA,*BtfC,ZO~ 

EN-N 

EN-F 

ZA-EN+1. 

-  'Zf»IH*EM*I«  -  _  _  . 

ZC-EN+EN-ENU 

Z0iCN*EM*l.*|N0 

EK«P$A|ZA)4P5*tZB)-PSAUC)-PSAfZDI 

bOurn 

END 

JlBFfC  EM. 

DOUBLE  PRECISION  FUNCTION  EKPtN»N#ENU) 

'  DOUBLE  PRECISION  FSB 

DOUBLE  PRECISION  EKPtEN,EN,ENU»ZA,ZB,ZC.ZO 
EN-N 
EN-P 

"  zC-en*en-Enu 

ZD-EN*EN*l.*ENU 
EKP-PSB(ZC)'PSBiZD) 

_ RETURN _  _ _ 

ENO~  . 

ftIBETC  6. 

_  - DOUBLE  PRECISION  FUNCTION  6(R.N.CWr 

DOUBLE  PRECISION  PSA 

OOUBLE  PRECISION  ENU»C,EN»EN'f ZA,ZB,ZC,ZD 

EN-N _  _ _ _ 

EN-N 

ZA-EN-ENU 
ZB-EN*1.*ENU 
ZC-EN-ENU 
-  ZD»EN*i.*ENU 

_  C  — PSA(ZAMPSA(ZBI*PSACZCI~PSACZD) 

RETURN  '  "  ~ 

ENO 
ftIBETC  M. 

OOUBLE  PRECISION  FUNCTION  HIN.N.ENUI 
DOUBLE  PRECISION  PSA 
OOUfLE  PRECISION  M.EN»EW,ZA»ZB,ZC»ZO»ENU 
EN-N 
EN-P 

ZA>EN»EK-ENU 

ZB«EN*EN*l.*ENU 

ZC-ER-ENU 

ZO-Eff+l.+ENU 

H~  —PSA  I Z Af *P  Sa  t  ZB I  ♦PSRTZC  r-RSA  CIO  T - 

RETURN 
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